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Abstract 

We analyze the stability and stabilizability properties of mixed retarded-neutral 
type systems when the neutral term is allowed to be singular. Considering an operator 
model of the system in a Hilbert space we are interesting in the critical case when 
there exists a sequence of eigenvalues with real parts approaching to zero. In this case 
the exponential stability is not possible and we are studying the strong asymptotic 
stability property. The behavior of spectra of mixed retarded-neutral type systems 
does not allow to apply directly neither methods of retarded system nor the approach 
of neutral type systems for analysis of stability. In this paper two technics are combined 
to get the conditions of asymptotic non-exponential stability: the existence of a Riesz 
basis of invariant finite-dimensional subspaces and the boundedness of the resolvent in 
some subspaces of a special decomposition of the state space. For unstable systems the 
technics introduced allow to analyze the concept of regular strong stabilizability for 
mixed retarded-neutral type systems. The present paper extends the results on stability 
obtained in [R. Rabah, G.M. Sklyar, A. V. Rezounenko, Stability analysis of neutral 
type systems in Hilbert space. J. of Differential Equations, 214(2005), No. 2, 391-428] 
and the results on stabilizability from [R. Rabah, G.M. Sklyar, A. V. Rezounenko, 
On strong regular stabilizability for linear neutral type systems. J. of Differential 
Equations, 245(2008), No. 3, 569-593]. Comparing with the mentioned papers, here 
we avoid a restrictive assumption of non-singularity of the main neutral term. 
Keywords. Retarded-neutral type systems, asymptotic non-exponential stability, sta- 
bilizability, infinite dimensional systems. 

Mathematical subject classification. 93C23, 34K06, 34K20, 34K40, 49K25. 



1 Introduction 

The interest in considering delay differential equations and corresponding infinite-dimen- 
sional dynamical systems is caused by a huge amount of applied problem which can be 
described by these equations. The stability theory of such type of systems was studied 
intensively (see e.g. [21 El [H])- Number of results was obtained for retarded systems, 
however an analysis of neutral type systems is much more complicated and these systems 
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are still studied not so deeply. We consider neutral type systems given by the following 
functional differential equation: 

j^[z{t)-Kzt] = Lzt + Buit), t>0, (1.1) 

where Zt : [—1, 0] C" is the history of z defined by zt{9) = z(t + 6). We assume the delay 
operator L : if ^([— 1, 0], C") C" to be linear and bounded, thus, it has the following form: 

Lf = ^° A2{d)ne) de + A,{d)f{e) d^, / g h\[-i, o], c"), (1.2) 

where ^2, A^ are n x n-matrices whose elements belong to L2([— 1, 0], C). We take the 
difference operator K in the form 

i^/ = A_i/(-l), (1.3) 

where A^i is a constant n x n-matrix. The form (11. 3p may be considered as a particular 
case of the operator K : C([-1,0],C") C" given by Kf = J^^ d^{9) f (9) , where /i(-) : 
[—1,0] — ^ i^nxn jg q£ bounded variation and continuous at zero. However, the considered 
model (II. 3p is sufficiently general. Its analysis is difficult enough and the results obtained 
are derived, in part, from the properties of the matrix A_i. 

The well-known approach, when studying systems of the form (II. ip . is to consider a 
corresponding infinite-dimensional model x = Ax, where A is the infinitesimal generator of 
a Co-semigroup. For systems (ll.ip -( |T73l) the resolvent of the operator A allows an explicit 
representation (see [Sni EI])- Such a representation is an effective tool for analyzing the 
exponential stability property since the last is equivalent to the uniform boundedness of the 
resolvent on the complex right half-plane. The resolvent boundedness approach is exhaustive 
when one considers stability of pure retarded type systems (A-i = 0) since such systems 
may be exponentially stable or unstable only. This fact is due to that exponential growth 
of the semigroup {e*'^}t>o is determinated by spectrum's location and there are only a finite 
number of eigenvalues of A in any half-plane {A : Re A > C}. 

For neutral-type systems (A_i ^ 0) in addition to the notion of exponential stabil- 
ity, which is characterized by the condition that the spectrum is bounded away from the 
imaginary axis (see [HI Theorem 6.1], [6]), one meets the notion of strong asymptotic non- 
exponential stability. This type of stability may happen in some critical case when the 
exponential stability is not possible (see e.g. [3]). Thus, strong stability cannot be described 
in terms of the resolvent boundedness. In [201 EI] for neutral type systems with a nonsingu- 
lar neutral term (det A_i ^ 0) this type of stability was precisely investigated for systems 
of the form (ll.ip -( fL3l) and some necessary and sufficient conditions of strong stability and 
instability had been proved. The proofs are based on such a powerful tool as existence of 
a Riesz basis of ^-invariant finite-dimensional subspaces of the state space and on further 
application of the results on strong stability in Banach spaces that had been originated in 

and later developed in [H dSl ESI El] and many others (see e.g. [2Z] for a review). 

In the case of neutral type systems with a singular neutral term (det A_i = and 
A-i 7^ 0), which we call mixed retarded-neutral type systems, the strong stability may also 
happen. However, the approach given in [201 EI] cannot be directly applied to such systems, 
since the existence of a Riesz basis of ^-invariant finite-dimensional subspaces of the whole 
state space cannot be guarantied. Moreover, mixed retarded-neutral type systems, in general, 
cannot be decomposed onto systems of pure neutral and pure retarded types. Therefore, the 
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analysis of strong stability for mixed retarded-neutral type systems poses a hard problem 
which requires bringing in essential ideas in addition. 

The method presented in this paper is based on a decomposition of the initial infinite- 
dimensional model X = Ax onto two systems xq = AqXq and xi = AiXi in such a way 
that the spectra of -^i satisfy: Re <7{Ao) < —e and —e < Re a{Ai) < for some e > 0. 
Generally speaking, the operators and Ai are not the models of delay systems (retarded or 
neutral), what, in particular, implies that the relation between their exponential growth and 
spectrum's location is unknown a priori. We prove the exponential stability of the operator 
^0 analyzing the boundedness of its resolvent. This direct analysis requires subtle estimates 
and the proof is technically complicated. For the analysis of the subsystem xi = AiXi we 
apply methods of strong stability introduced in [201121]. Finally, the introduced approach 
allows us to prove for mixed retarded-neutral type systems the results on strong stability 
formulated in [2T] . 

Besides, for control systems the proposed approach allows to analyze the notion of 
regular asymptotic stabilizability [22] which is closely related to the strong stability notion. 
The technic of the regular asymptotic stabilizability were introduced in |22j and the sufficient 
condition for the system fll.ll) - fll.3l) to be stabilizable had been proved in the case det A_i ^ 
0. In the present paper, using the same framework as for stability, we show that these results 
hold for mixed retarded-neutral type systems also. 

The general framework which we use is the theory of Co-semigroups of linear bounded 
operators (see e.g. [2Z!)- In order to precise the main contribution of our paper let us first 
give the operator model of the system fll.ip - fll.3l) . We use the model introduced by Burns 
et al. [1] in a Hilbert state space. The state operator is given by 

Ax(t) -a( y^^^ ] - ( ^2(^^)it(^^) d^ + A,{o)zm de \ 

with the domain 

V{A) = {(y, z{-)f : z e H\-l, 0; C"), y = ;.(0) - A_iz(-1)} C M2, (1.5) 

where M2 = C" x L2(-1,0;C") is the state space. The operator A is the infinitesimal gen- 
erator of a Cq- semigroup. Studying stability problem, we consider the model 

X = Ax, x{t) = ( ] , (1.6) 



corresponding to the equation fll.ll) - fll.3p with the control u = 0, i.e. to the equation 

A2{9)z{t + 9) d9 + J A3{9)z(t + 9) d9, t>0. (1.7) 

The solutions of (11. 7p and (11.60 are related as Zt{9) = z{t + 9), 9 E [—1, 0]. 
For the control the equation (ll.ip -( |L3l) which we rewrite as 

z{t)=A_iz{t-l)+ I A2i9)z{t + 9)d9+ I A3{9)z{t + 9) d9 + Bu, (1.8) 

we consider the model 



-1 j-i 



x = Ax + Bu, (1.9) 
where the operator B : ^ M2 is defined by n x p-matrix B as follows: Bu'= {Bu, 0)^ . 
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The operator A given by (ll.4p possesses only discrete spectrum cr{A), and, moreover, the 
growth of the semigroup {e*"^}t>o is determinated by spectrum's location. Namely, denoting 
by Ug = sup{ReA : A G o-{A)} and by Uq = mi{u : ||e-^*x|| < Me'^*||x||}, we have the 
relation uq = ujg (see e.g. [6j|). 

For stability problem, the last fact implies that the semigroup {e*"^}t>o is exponentially 
stable if and only if the spectrum of A satisfies cu^ < 0. However, this type of stability is not 
the only possible one for systems of the form (11.71) (e.g. the same situation can also happen 
for some hyperbolic partial differential equation). Namely, if tUs = (and A^i ^ 0), then 
there exists a sequence of eigenvalues with real parts approaching to zero and imaginary part 
tending to infinity. In this critical case the exponential stability is not possible: He*"^! -/^ 
when t — > oo, but asymptotic non-exponential stability may occur: lim e*-^x = for all 

t— >+oo 

X G M2. For systems (11. 6p . satisfying the assumption dety4_i 7^ 0, the problem of strong 
stability was analyzed in ^20 ^T\. The main result on stability obtained there may be 
formulated as follows. 

Theorem 1.1 ([211 R- Rabah, G.M. Sklyar, A.V. Rezounenko]). Consider the system 1^1. 6\) 
such that det A-i 7^ 0. Let us put a\ = cr(A_i) fl {fi : = 1}. Assume that o-{A) C {A : 
ReA < 0} (necessary condition). The following three mutually exclusive possibilities hold 
true: 

(i) (Ji consists of simple eigenvalues only, i.e. an one- dimensional eigenspace corresponds 
to each eigenvalue and there are no root vectors. Then system 1^1. 6\} is asymptotically stable. 

(ii) The matrix A^i has a Jordan block, corresponding to an eigenvalue /i G (Xi. Then 
lll.6\) is unstable. 

(ii) There are no Jordan blocks, corresponding to eigenvalues from ai, but there exists 
an eigenvalue G o"i whose eigenspace is at least two dimensional. In this case system lil.6\) 
can be either stable or unstable. Moreover, there exist two systems with the same spectrum, 
such that one of them is stable while the other one is unstable. 

Let us discuss the importance of the assumption det A^i 7^ 0. The proof of Theorem ll.il 
given in is based on the following facts. Firstly, if det A_i 7^ 0, then the spectrum of 
A is located in a vertical strip di < Re (t{A) < d2. Namely, in [Sni El] it was shown that 
(j{A) = {In \firn\ + i(arg/im + 2Trk) + o{l/k) : fim G a{A^i), k G Z}. From the last it also 
follows the necessary condition for the system to be asymptotically stable: a{A^i) C {/i : 

l/^l < !}• 

Secondly, such location of the spectrum had allowed to prove the existence of a Riesz 
basis of generalized eigenvectors for the operator A = A corresponding to the case ^2(6*) = 
A^IO) = 0. For a general operator A the generalized eigenvectors may not constitute a basis 
of the state space (see an example in [20] and some general conditions in [28j). However, 
in [2ni [21] it was proved the existence of a Riesz basis of ^-invariant finite-dimensional 
subspaces of the space M2 (see also [22] )• Such a basis is a powerful tool that had been 
apphed for the analysis of strong stabihty. 

If we allow the matrix A^i to be singular, then the described above location of the 
spectrum does not hold anymore. Generally speaking, in this case for any a G M there exists 
an infinite number of eigenvalues which are situated on the left of the vertical line ReA = a. 
Thus, the existence of a Riesz basis of ^-invariant finite-dimensional subspaces for the whole 
space M2 cannot be guarantied. As a consequence, the proof of the item (i) given in [21], 
which is essentially based on the Riesz basis technic, is no longer satisfactory and one needs 
another way of the analysis of stability. 

However, it can be asserted that nonzero fim ^ cr(^-i) define the spectral set {In \ firn \ + 
i(arg/im + 27rA;) +o(l/A;) : fim ^ o"(^-i),/^m 7^ 0, /c G Z} C o-{A) which belongs to a 
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vertical strip. In particular, this can be asserted for /x^ G cti. The fact that Theorem 11.11 
is formulated in terms of ai and the last remark give us the idea to decompose the initial 
system (11. 6p into two systems 

I Xi = AiXi 

in such way that <7{Ao) = cri-A.) fl {A : — oo < ReA < —e} and <7{Ai) = <7{A) fl {A : —e< 
ReA < ujs = 0}, for some e > 0. 

To obtain the representation (II.IUI) we construct a special spectral decomposition of 
the state space: M2 = ® Ml, where Mg, Ml are ^-invariant subspaces. We define 
Aq = A\mo and Ai = A\m^- 

The spectrum of the system xi = AiXi is such that the corresponding eigenvectors form 
a Riesz basis of the subspace Ml. The strong stability of the semigroup {e*'^|jv/i}t>o is being 
proved using the methods of pTj . 

The semigroup {e*'^|jv/o}t>o is exponentially stable. We prove this fact using the equiva- 
lent condition consisting in the uniform boundedness of the resolvent R{\, A)\mo on the set 
{A : ReA > 0}. Thus, we prove that the initial system x = Ax is asymptotically stable. The 
mentioned scheme requires complicated technics. 

To complete the stability analysis we revisit the example showing the item (iii) with 
a simpler formulation than in [21], where it was given using the Riesz basis technic. The 
analysis of the spectrum being carried out in our example is essentially based on the deep 
results on transcendental equations obtained by L. Pontryagin ^7\. We notice also that the 
proof of the item (ii) given in [2T] does not involve the Riesz basis technic and, thus, it 
remains the same for the case det A_i = . 

The technics of the direct spectral decompositions and the resolvent boundedness pre- 
sented above allow us to extend the results on the stabilizability problem given in [22] for 
the case of singular matrix 

The general problem of stabilizability of control system is to find a feedback u = Tx 
such that the closed-loop system 

x = (^ + BT)x 

is asymptotically stable in some sense. For the system (II. 8p the result of exponential stabi- 
lizability may be derived from those obtained for some particular cases (see e.g. [71 [TSl [T6] ). 
The needed feedback for our system is of the form 

/o /-O 
F2(^)i(t + d^ + y F^{e)z{t\Q)&e. (1.11) 



Our purpose is to obtain, as in [22], the condition of asymptotic non-exponential stabiliz- 
ability of the system (11.81) with the regular feedback 



/o /-O 
F2i9)z{t + 9) d9 + J F3{9)z{t + 9) d9, (1.12) 

where F2{-), F^l-) G L2(— 1, 0; C"^^). The motivation is that this kind of feedback is relatively 
bounded with respect to the state operator A and does not change the domain of A: T>{A) = 
T>[A+BT). The natural necessary condition regular stabilizability is (J (A_i) C {/i : < 1} 
because A_i is not modified by the feedback. Under the same restrictive condition det A_i 7^ 
in [22j was obtained the following result on stabilizability. 
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Theorem 1.2 ([221 R- Rabah, G.M. Sklyar, A.V. Rezounenko] ) . Let the system l\1.8]) verifies 
the following assumptions: 

(1) All the eigenvalues of the matrix satisfy < 1. 

(2) All the eigenvalues /x G (Xi are simple. 

Then the system is regularly asymptotically stabilizable if 

(3) rank(A(A), B) = n for all A : ReA > 0. 

(4) rank(/i/ — B) = n for all /i G ai. 

The proof of this theorem given in [22] uses the existence of the Riesz basis of the whole 
state space M2 and, thus, it requires the assumption detA_i 7^ 0. To avoid this assump- 
tion, we construct and prove another spectral decomposition which takes into account the 
unstable part of the system. By means of this decomposition we separate a subsystem which 
is generated by the part of the spectrum corresponding to the zero eigenvalues, i.e. the 
singularities of the matrix A_i. Proving the resolvent boundedness, we show the exponen- 
tial stability of this subsystem. The main "critical" part of the system is in ^-invariant 
subspaces, where we apply the same methods that were given in [22], namely, the theorem 
on infinite pole assignment, introduced there, and a classical pole assignment result in finite 
dimensional spaces. 

The paper is organized as follows. In Sections 2 we recall the results on the spectrum, 
eigenvectors and the resolvent of the operator A obtained in [^ [22] • Besides we prove 
some properties of eigenvectors. In Section 3 we construct and prove two direct spectral 
decomposition of the state space. One of them is used to prove the main result on stability 
and another one for the proof the result on stabilizability. Section 4 is devoted to the proof 
of the uniform boundedness of the restriction of the resolvent on some invariant subspaces. 
Finally, in Section 5 and Section 6 we give the formulation and the proof of our main results 
on stability and stabilizability. Besides, in Section 5 we give an explicit example of two 
systems having the same spectrum in the open left half-plane but one of these systems is 
asymptotically stable while the other one is unstable. 

2 Preliminaries 

In this section we recall several results on the location of the spectrum of the operator A, 
on the explicit form of its resolvent and on the form of eigenvectors of A and ^4*. We prove 
some properties of eigenvectors of A and A* . 

2.1 The resolvent and the spectrum 

The results given in this subsection have been presented and proved in [^ [211 [22]- Some 
formulations of the propositions are adapted for the case dety4_i = 0. 

Proposition 2.1 ( [2T| Proposition 1]). The resolvent of the operator A has the following 
form: 



\^^'> ^ \ e^(^-^)e(s)ds + e^^A-i(A)D(2,^,A) / 

where z G C",^(-) G L2(— 1,0;C"); A(A) is the matrix function defined by 

A(A) = A^(A) = -AJ + Ae-M_i + A / e^'A2{s)ds+ e^'A3{s)ds, (2.14) 




e"M_i e-^'^{s)ds + (/ - e-M_i) A-i 




-1 J-i 
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and D{z,C,, A) is the following vector-function acting to C"; 

r-O 



D{z, A) = z + Ae-M_i / e-^'^iO) dO 







-1 



-As 



as) ds 



de. 



(2.15) 



From f l2.13p one may see that the resolvent does not exist in the points of singularity of 
the matrix A (A), i.e. the equation det A(A) = defines the eigenvalues of the operator A. 
Now let us characterize the spectrum of A more precisely. 

We denote by /ii, . . . , /i£ the set of distinct eigenvalues of the matrix A_i and hj pi, . . . ,pi 
their multiplicities. We recall the notation ai = cr{A_i) fl {/i : |/i| = 1} and assume that 
0"! = {/^i, • • • , yU^i}, ^1 < ^- We notice that one of the eigenvalues yU£i+i, ■ ■ ■ , fii may be zero. 

Further, studying stability and stabilizability problems, we consider mainly the situa- 
tions when the eigenvalues from ai are simple. This gives us a motivation to assume below 
(if the opposite is not mentioned) that pi = . . . = = 1. Besides, without loss of generality, 
we assume that the matrix A_i is in the following Jordan form: 



/ /^l ■ 


. 







. 







. 


. 


. 


Jii+i ■ 


. 


V . 


. 





. Ji J 



(2.16) 



where J^^+i, ■ ■ ■ , Ji are Jordan blocks corresponding to the eigenvalues /i^^+i, . . . , 

Let us denote by A the state operator in the case when A2{0) = A3 (6*) = 0. It is not 
difficult to see that the spectrum of A has the following structure 

a{A) = {A^ = In + i(arg/im + 2nk) : m = 1, . . 



/i„^0,A;eZ}U{0}. 



We denote by L'^ir^^^) circles centered at A,^ with radii r'^^\ 



Proposition 2.2. Let ai = {/ii, . . . , /i^^} consists of simple eigenvalues only. There exists 
A^i e N such that the total multiplicity of the roots of the equation det A (A) = 0, contained 
in the circles L^(r(^)), equals = 1 for all m = 1, . . . , £1 and k : \k\ > Ni, and the radii 
r*^^^ satisfy the relation J2{r^''^Y < 00. 

This proposition is a particular case of [221 Theorem 4] which have been formulated and 
proved there under the assumption det A^i 7^ 0. The proof for the case det A_i = remains 
the same. 

Notation 2.1. We denote the eigenvalues of A mentioned in Proposition by X^, m = 
l,...,£i, \k\>N,. 



Remark 2.3. Proposition \2.2\ is formulated for m = 1, . . . however, it also holds for all 
those indices m = 1, . . . ,i which correspond to nonzero eigenvalues fim ^ cr(A-i)- 

Remark 2.4. In the case detA_i 7^ the spectrum of A belongs to a vertical strip which 
is bounded from the left and from the right. However, in the case det A_i = 0, in addition 
to the eigenvalues mentioned in Proposition \2.S\ the operator A may also possess an infinite 
sequence of eigenvalues with real parts tending to —00. 
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Similar results hold for the operator A*. The spectra of A and A* are related as 
cf^A*) = a{A). Eigenvalues of A* are the roots of the equation det A*(A) = 0, where 

/O pO 
e^'Al{s)ds + J e^'Al{s)ds, (2.17) 

and the relation (A(A))* = A*(A) holds. The eigenvalues A^, m = 1, . . . \k\ > Ni may 
be described as in Proposition 12.21 



2.2 Eigenvectors of A and A* 

First we give the explicit form of eigenvectors which has been proved in 

Proposition 2.5 ([^ Theorem 2], [221 Theorem 7]). Eigenvectors (f: {A — XI)<f = and 
tp: {A* - = 0, A e (t{A) are of the form: 



^ = <^(A) 



(/ - e'^A^i)x 



^ = V(A) = 



Ae-^^ - A;{e) + e-^^ / e^^(v4*(s) + XA;{s))ds 



(2.18) 
(2.19) 



where x = x(A) G KerA(A), y = y{\) G KerA*(A). 

Below we give several properties of the sets of eigenvectors and we begin with the 
calculation of the scalar product between eigenvectors of A and A*. 

Lemma 2.6. Let Aq, Ai G (y{A) and (p = f{Xo), "ip = V'(Ai) are corresponding eigenvectors: 
{A — XqI)^ = 0, {A* — = 0. The scalar product {(f, iP)m2 ^Quals to the following value: 



Ah 



0, Ao ^ Ai 

-(A'(Ao)a;,y)c", Aq = Ai 



(2.20) 



where A'(A) = ^A(A) and x = x(Ao), y = l/(Ai) are defined by Ii2.18\) and ^2.1fA) . 

Proof. First let Aq 7^ Ai and we compute directly the scalar product (v?, ip) M2 using the 
representations f l2.18p and (12.191) : 



M2 



((/ - e-^M_i)x, y)c" + / (e^o^x, \ie-^^'y)c- d9 - J {e^^'x, A*(^)i/)c" d^ 

-1 -1 

+ / /e^o^x,e"^^/e^"(A*(s) + A7A*(s)) ds ■?/) d^ 

-1 \ / C" 

((/ - e-^M_i)x, y) + (^I^ Aie(^»-^i)^ dO ■ x,y'^ - / 6^0^2(0) d^ ■ x, 

/ e \ 

+ / / e(^«-^i)^ / e^''{A^{s) + AiA2(s)) ds d^ ■ x, y) 



(r(Ao,Ai)x, y) 



where 



r(Ao,Ai) = J-e-^M_i + Ai / e(^"-^i)^ d^- / e^oM2(0) d^ 

-1 -1 

e 
+ J e{\o-\i)e J e^is^Asis) + \iA2{s)) ds d9. 



(2.21) 



(2.22) 



The last term of r(Ao, Ai), which is the integral over the domain — 1 < < s < 0, we rewrite 
using the identity J^^ G{s, 6) ds d6' = — J°-^ J^-^ G(s, 0) AO ds which holds for any function 
G{s,6). Taking into account the relation J^^ei^o->^i)s _ _ gAi-Ao^^ obtain 

e s 

J e(^o-^i)''/e^i^(A3(s) + AiA2(s)) ds d^ = - / e^^'{As{s) + AiA2(s)) / e^^"-^^)^ d^ ds 

-1 -1 -1 



gAi-Ao J e^i«(A3(s) + AiA2(s)) ds - / e^«^(A3(s) + AiA2(s)) ds . 
-1 -1 



Aq — Ai 



Finally, we have 
r(Ao,Ai) 



[(Ao - Ai)/ - (Ao - Ai)e-^M_i + Ai(l - e^^-^o)/ 



Aq — Ai 



-(Ao - Ai) / e^^'A^ie) de- J e^«^(A3(s) + AiA2(s)) ds 



-1 



"1 







+gAi-Ao J e^i^(A3(s) + AiA2(s)) ds 



1 



Aq — Ai 



Ao/ - Aoe~^M„i - / e^o^(A3(s) + AoA2(s)) ds 

-1 



-Aie^i-^oj + Aie-^M_i + e^^-^" / e^^'{A3{s) + AiA2(s)) ds 

= A;rAr[-^(^o)+e^-^oA(AO]. " 

Taking into account that x G KerA(Ao), y G KerA*(Ai) and (A(Ai))* = A*(Ai), we 
conclude that 

„Ai — Ao 



M2 



1 



Ao — Ai 



[-A(Ao)+e^-^»A(Ai)] x,y 



-n Aq — Ai 



(x, A*(Ai)i/)cn = 0. 



(2.23) 



Let us now consider the case Aq = Ai. From (12.211) . (I2.22p we have: 

{(P,^P)m2 = (r(Ao)x,?/)c", 

where 



r(Ao) = / - e^^M_i + Ao/ - / e^«^A2(^) d^ + / / e^"^(A3(s) + AoA2(s)) ds d^^. (2.24) 

J-i J~i Jo 

The last term of r(Ao), which is the integral over the domain — 1 < < s < 0, we 
rewrite using the identity J^^ G{s, 6) ds d6' = — j^^ J^^ G(s, 6) d6 ds. Thus, we obtain: 

r(Ao) = /-e-^oA_i + Ao/- /e^«M2(^^)de- /e^o^(A3(s) + AoA2(s)) / d^ds 

-1 -1 -1 

= (^/ - e-^M_i - / e^«^(sA3(s) + sAoA2(s) + A2(s)) ds^ 

+ (^\oI - / e^«^(A3(s) + AoA2(s)) ds^ 

= -A'(Ao)-A(Ao). 
Taking into account the relation x G KerA(Ao), we conclude that 

{ip,^)M2 = -{A'i\o)x,y)c^. (2.25) 
The last completes the proof of the lemma. □ 
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For v^(A^) and ip{X^) we will use the notation ip^ and ip'^ respectively. Besides, we use 
x'^ and instead of a;(A^) and ?/(A^). 

Lemma 2.7. Let a\ = {yUi, . . . , /i^^} consists of simple eigenvalues only. The eigenvectors 
ip^, m = 1, . . . , Ci, k : \k\ > Ni constitute a Riesz basis of the closure of their linear span. 
The same holds for eigenvectors ipm' = 1, • • • , ^i? k : \k\ > Ni. 



A more general formulation of this proposition have been given in [20i Theorem 7, The- 
orem 15] under the assumption detA^i 7^ 0. We give a sketch of the proof in our case. 

The families of functions {e^"^^}k(^z form an orthogonal basis of the space L2([— 1, 0], C) 
for each m = 1, . . . , £1, where = i(arg /i^ + 2'7ik) are eigenvalues of the operator A. Thus, 
the functions {G^'^^}\k\>N, N eN form a basis of the closure of their linear span. 

Since we have chosen the matrix A_i in the form 02.161] and due to (12.181) . the eigen- 
vectors if'L of A are of the form i^jL = ( „ ) , = (0, . . . , 1, . . . , 0)"^. Therefore, the 

family : m = 1, . . . , £1 : > A^^i} is a basis of the closure of its linear span. 

/ (/ - e~^™A_i)x^ \ 

The eigenvectors (pz^ = \k a i of ^ are quadratically close to (pz^. To 

V e-'-'^x^ J 

prove this fact we should argue similar to Theorem 15 given in [20] (see also ^Q\)- Thus, 
eigenvectors ip^, m = 1, . . . , ii, k : \k\ > Ni constitute a Riesz basis of the closure of their 
linear span. 

From Lemma 12.61 and Lemma 12.71 we conclude the following. 

Corollary 2.8. The sequences and ip^, m = l,...,£i, k : \k\ > Ni are biorthogonal 
after normalization and ^^)m2 = -(A'(A^)x^, ?/^)c". 

The following relation will be essentially used in the analysis of the boundedness of the 
resolvent in Section 4. 

Lemma 2.9. Let ip = iIj{Xq), Aq G cr(^) be an eigenvector of the operator A* and let 
g = (z,^(-)) G M2 be orthogonal to ip: g-Lip. Then the following relation holds: 

D(^,e,Ao) GlmA(Ao), (2.26) 

where D{z,$,,\) is defined by ^2.1\) . 

Proof. We show the relation D{z, ^, Ao)-LKerA*(Ao) which is equivalent to (I2.26p . The eigen- 
vector ip is of the form (12.191) : 



y 



Aoe-^«^ - Al{e) + e-^o^ J^e^'>'Al{s) ds + Aoe"^«^ /o e^°M*(s) ds 



y 



where y = y{Xo) G KerA*(Ao). For any g = {z,^{-)), which is orthogonal to ip, we obtain: 

= {g,i^)M, = {z, + / (mAoe-^"'y) d9-J im, ^2(^)z/)c" 



+ / UiP). e-^"'^ / e^°^(A*(s) + AoA*(s)) ds ■ i/\ d^ 

-1 \ / C" 

\-l / c" \-i 





(2.27) 



-1 



m d9, y 
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Using the identity J^^ G{s, 9) ds dO = — J^^ J^-^ G{s, 9) dO ds which holds for any function 
G(s, 6*), we rewrite the the last term of (12.271) . and, finally, we obtain the relation: 



(2.28) 



\ -1 -1 

- J e^o^ [A,{s) + XoA2{s)] J e-^«^^(0) d^^ ds, y 
-1 -1 

Since y G KerA*(Ao), then for any x G C" we have: 

= (a;, A*(Ao)?/)c" = (A(Ao)a;, y)c^. 

Therefore, for any 6 the relation (e~'*'°^A(Ao)^(^^), y)c" = holds, and, integrating it by 6 
from —1 to 0, we obtain: 

= (| / e-^«^A(Ao)e(^) dO, = / Aoe~^o^e(^) dO + / Aoe-^«M_ie(^) d^ 

+ / e^o^ [^3(5) + AoA2(s)] ds / e-^«^e(^) ^6, y 

—1 —1 / 

(2.29) 

Let us sum up the left-hand sides and the right-hand sides of the relations (12.291) and (12.281) . 
In the obtained relation the term J^^ Aoe~'^''^^(6') dO is cancelled. The last terms of (12.291) and 
(I2.28P we sum up according to the identity — /^^^ G{s, 6) dO ds + G{s, 0) d6 ds = 

- /_i So ^) dO ds = - j\ G{e, s) ds de which holds true for any function G{s, 6). 
Finally, we obtain: 

/ 

= lz + \oe-^^A_ije-^^^i{e)de-jA2{e)i{e)de 
\ ~i -1 





- / e^o^ [A,{e) + A0A2 



-1 



/ e-^«^^(s) ds 
.0 



d^, y 



Since y G KerA*(Ao), we conclude that i5(2;, ^, Aq) J-KerA*(Ao), what completes the 
proof of the lemma. □ 

Remark 2.10. We emphasize the fact that det A(Ao) = and, therefore, the matrix A^^{Xq) 
does not exist. However, the proved relation D{z,^,Xq) G ImA(Ao) means that there exists 
the inverse image of the vector D{z,C,, Aq) with respect to the matrix A(Ao). 



3 Spectral decompositions of the state space 

We recall that we consider the operator A in the case when all eigenvalues from (Xi C cr(y4_i) 
are simple. In this section we construct construct two spectral decompositions of the state 
space M2. Assuming that (j{A) C {A : ReA < 0}, in the first subsection we construct 
a decomposition which we further us in Section 5 for the stability analysis. In the second 
subsection we assume only < 1 for all fi G a{A^i) (i.e. a part of the spectrum of A may 
belongs to the closed right half-plane) and construct a decomposition needed in Section 6 
for the stabilizability analysis. The structures of these decomposition are very similar. In 
the third subsection we prove some technical results used in the proofs of validity of the 
decompositions. 
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3.1 Spectral decomposition for the stability problem 



For the stability analysis our aim is to divide the system onto exponentially stable part and 
strongly asymptotically stable part. To do this we construct a decomposition of the state 
space M2 onto the direct sum of two ^-invariant subspaces and prove its validity. 
We divide the spectrum of A onto two parts. For some N > Ni we define 

Ai = Ai(iV) = {A^ G a{A), m = 1, . . . Julkl > N}, (3.30) 

and represent the spectrum as follows: 

a{A) = AqU Ai. 

Remark 3.1. The set Ai is determined by N ^ N. For any small e > there exists big 
enough N such that Ai belongs to the vertical strip {A : — e < ReA < 0}. 

The following figure illustrates our idea: 




© 




© 





© 



© 

® 



© 



Figure 1. 




Re 



By crosses we denote A^, eigenvalues of A and by points we denote eigenvalues of the 
operator A. 

We introduce two subspaces of M2: 

= M^{N) = CI Lm{(p : {A - = 0, A G AJ, 



(3.31) 

'2 ~ ^"2 V — -^"nV" ■ V-^ — /v-i ;ity — w, /V c iilj. (3.32) 

Obviously, is ^-invariant and is ^*-invariant. We introduce = M2{N) which 
satisfies 



M} = M^{N) = CI Lin{V^ : (A* - A/)V^ = 0, A G AJ. 



M2 = M^®M^. 
Due to the construction, M2 is an ^-invariant subspace. 



(3.33) 
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Remark 3.2. We recall that due to Lemma \2. 7| eigenvectors {v^^} of A, corresponding 
to AJ^ G Ai, form a Riesz basis of the closure of their linear span. The same holds for 
eigenvectors {V'^} of A* , corresponding to G Ai. 

The main result of this subsection is the following theorem. 

Theorem 3.3 (on direct decomposition). Let ai = {/xi, . . . , fie^} consists of simple eigenval- 
ues only. For any N > Ni the subset Ai = Ai(A^) C cr(^) given by hS. 30\) and the subspaces 
M2, Ml, Ml, given by ^3. 31]) . ^3. 32^ and ^3. 33\) define the direct decomposition of the space: 

M2 = Ml © M° (3.34) 

where the subspaces Ml, M2 are A-invariant. 

Proof. To prove (13.341) we show that any element ^ G M2 allows the following representation: 
^ = ^o + EE^-^- ^0 e M°, y.^GM2\ 5^ 5^ |4P< 00. (3.35) 

m=l |fc|>Af m=l |fc|>Ar 

As we have noticed above, the eigenvectors {99^ : {A — X'^I)ip^ = 0, A^ G Ai} given 
by (12.181) form a Riesz basis of the closure of their linear span. Besides the eigenvectors 
{=7/'^ : {A* — X^I)'iIj^ = 0, A^ G Ai} form a Riesz basis of the closure of their linear span, 

where eigenvectors are given by (I2.19p . We use the notation '(p^ = =V'm- 

Using the representations (12.181) and (I2.19p . we choose eigenvectors with ||a;^||c" = 1 
and lll/mllc" = 1- Due to Lemma [331 given in Subsection 13.31 there exists C > such that 
llv^mlU/^ < C and H^^IU/a < C for all m = 1, . . \k\ > N. 

Applying the decomposition (I3.33P to vectors ip^, we obtain 

1=1 \j\>N 

Since (v?^, ipj) = for (m, k) 7^ (i, j) (Corollary 12. 81) . the last representation may be rewritten 
as follows: 

vi = 7t + a'Jt, 7™eM°, (3.36) 
moreover, due to (I2.20p we have the relation 

^ {vL^'Jm, ^ i^i^'m.i^'jM, ^ -^(AXA^)x^,y^)cn 

From (I3.36P and (13.370 it also follows that 



h^ll < l|^^l| + |aj;,|l|z?^|| <C+./|-^(A'(A^)xLl/^)|. (3.38) 

y m 

Using the decomposition (I3.33P and the relation (I3.36p . we represent each vector ^ G M2 
as follows: 



e = eo+ E E hliii = ^,- E E ^7i+ E E Wm 

»n=l|fe|>Af m=l|A:|>Ar ™ m=l\k\>N 

= eo + E E cl^l. 

m=l |fc|>Ar 
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(3.39) 



where foe MO, ± E l&^P < oo, = - E E 1^7^ e M", 4 

m=l|fc|>Af m=l\k\>N 



To prove the vahdity of the decomposition fl3.39p it is enough to show that 

and ||7mll ^ C2 for some < Ci,C2 < +00. Taking into account (13 .3 71) and fl3.38p . the last 
means to give the estimate 



< Ci < 



A 



(3.40) 



This estimate is proved by Lemma 13.111 given in Subsection 13.31 

Thus, the representation (13.391) holds for any ^ G M2, what completes the proof of the 
theorem. □ 



3.2 Spectral decomposition for the stabilizability problem 

We recall that for the stabilizability problem we assume only that < 1 for all /i G a{A_i). 
In this case an infinite number of eigenvalues may belong to the right-half plane. On the 
other hand, only a finite number of eigenvalues may be located on the right of a vertical 
line ReA = e for any e > 0. For the analysis of stabilizability it is convenient to construct a 
decomposition of the state space onto three .4-invariant subspaces. 
We divide the spectrum of A onto three parts: 

a{A) = Ao U Ai U A2, (3.41) 

where the subsets Aq, Ai, A2 are constructedjjy the following procedure. 

Let A^o be such that A^ e L'^{r), r < 1|A^ - A^|, (m, k) ^ for all k > Nq and for 
all m such that Hm 7^ 0. First, we construct an auxiliary division 

a{A) = xi U xo, Xi = {A^ e ct{A) : \k\ > No,m = 1, . . . , i, fXm ^ 0}. 

Due to the construction, any vertical strip St(5i,(52) = {A : 61 < ReA < ^2} contains 
only a finite number of eigenvalues from xo- We also recall that Us = sup{ReA : A G 
a{A)} < +00. 

If (Ti 7^ then for any r > the strip St(— r, r) contains an infinite number of eigenvalues 
from xi and, as we have noticed above, only a finite number of eigenvalues from xo- Let us 
fix some r > and consider the value 



e = min |ReA 

AeSt(-r,r)nxo 



If £ > 0, then the vertical strip St(— £,£:) does not contain eigenvalues from xo and 
contains an infinite number of eigenvalues from xi- Moreover, the strip St(e,r) contains 
only a finite number of eigenvalues from xi- Thus, the strip St(£:,ti;s) contains a finite 
number of eigenvalues of the operator A and, therefore, we conclude that these eigenvalues 
are located in a rectangle {X : e < ReA < ujq, |ImA| < M} for some M > 0. Finally, we put 

Ag = a(A) n {A : ReA < 

Ai = (tU) nSt(-e,e), (3.42) 
A2 = a{A) n St{e,uJs). 

We illustrate the mentioned construction by the following figure: 
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operator A. 



We notice that the relation e = means that there exists eigenvalues with zero real 
part. In these case we calculate min |ReA| without taking these eigenvalues into 

AeSt(-r,r)nxo 

consideration and after constructing (13.421) we add these eigenvalues to A2. 

The obtained sets of eigenvalues may be described as follows: Aq belongs to the left 
half-plane and is separated from the imaginary axis; Ai consists of infinite number of simple 
eigenvalues which may be as stable as unstable, the corresponding eigenvectors form a Riesz 
basis of the closure of their linear span; A2 consists of finite number of unstable eigenvalues. 

Passing over to the construction of invariant subspaces, let us denote the elements of 
the finite set A2 as Aj, i = 1, . . . , r. We denote the corresponding generalized eigenvectors 
by (fij: {A — Xiiyifij = 0, j = 0, . . . , Sj — 1. As before, the eigenvalues from Ai we denote 
as and the corresponding eigenvectors we denote as ip'^, m = 1, . . . \k\ > N. 

We introduce the following two infinite-dimensional subspaces of eigenvectors: 

M', = CI Lin{¥,^ : {A - Xtl)^t = 0, Xt e AJ, 

Ml = CI Lin{^^ : (A* - A^/)^^ = 0, A^ G AJ, (3.43) 
two finite-dimensional subspaces of eigenvectors and root-vectors: 

= Uia{^ij : {A - Xiiy^^j = 0, Ai G A2, J = 0, . . . , Si - 1}, 

= Un{iJi,, : {A* - XiyiJij = 0, Ai G A2, J = 0, . . . , Si - 1} (3.44) 
and the subspace M2 , which satisfies 

M2 = (M^ © M|)©M°. (3.45) 
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Thus, we have constructed three ^-invariant subspaces: Mg, and M|. The main resuh 
of this subsection is the following theorem. 

Theorem 3.4. Let cxi = {/ii, . . . , fie^} consists of simple eigenvalues only. For any N > Ni 
the decomposition of the spectrum ^^-4^ and the subspaces given by ^EIM), ^EIM) and 
define the direct decomposition of the space M2: 

M2 = M° © © Ml 

where the subspaces M2, M^, M| are A-invariant. 

Proof. The proof of this proposition is very similar to the proof of Theorem 13.31 We prove 
that any element ^ G M2 allows the representation: 

ll rsi — l £1 

^ = ^0 + E E ^'rn^'rn + E E '^'^^^'^^ ^ ^2 , E E l^-l' < (^-^S) 

m=l \k\>N i=l j=0 m=l \k\>N 

The eigenvectors {y?^ : € Ai} form a Riesz basis of the closure of their linear span. 
The finite set of the generalized eigenvectors {(fiij '■ Aj G A2} is also a basis of their linear 
span. These vectors are linearly independent from eigenvectors {f^}- Thus eigenvectors 
{^p^} U {ipij} form a basis of the closure of their linear span. Arguing the same way, we 
conclude that the set {V'^} U {i'lj} is also a basis of the closure of its linear span, where 
= TT^'m ipij = T'^id- Morcovcr, without loss of generality, we may assume that the 

sets {v^m} U {i'm} U {i^ij} biorthogonal after the normalization. 

Using the representations fl2.18p and fl2.19p . we choose eigenvectors with ||a;^||c" = 1 and 
ll?/mllc" = 1- Due to Lemma 1331 there exists C > such that Hv^mlU/a < C and ||'?/'^||m2 < C 
for allm = \k\ >N. 

Applying the decomposition fl3.33p to vectors ip'^, we obtain 



£1 r Si-1 



2 • 



= 7™ + E E + E E ^M'V'^.i, e m; 

i=l |j|>Ar i=l j=0 



Since the sets {(f} and are biorthogonal, the last representation can be rewritten as 
follows: 

^t = lt + 7™ e Ml (3.47) 
and, moreover, due to Lemma [2.61 we have 

Besides, arguing the same we obtain: 

= 7i,i + (^iJ^iJiJ li,3 ^ M^ji = 0, . . . , Si - 1. 

From (13:471) and ([MHD it also follows that 



17^11 < ll^^ll + la^lll^^ll <C+\l\j{A'{Xi)xt,yi)\. (3.49) 
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Using the decomposition fl3.33p and the relation fl3.47p . we represent each vector C, € M2 
as follows: 

ii r Si~l 

^ = ^0 + E E + E E ^mV^m- 

m=l \k\>N 1=1 j=0 



in 



m=l \k\>N 



m 



^?,.7 . .... . - . .. . ^ ^7, 



i=l j=o *J 



m=l |fc|>Ar " ™ 

^0 + E E ^^'^^ + E E ^^j'^^'i' 



1 7=0 



m=l |fc|>Af 



i=l j=0 



where e M^, £ E l&^nP < 00, = ^o- E E S^7^-E E S^7., e M^, c„ - ^ 



r Si — 1 



(3.50) 



m=l |fc|>Af m=l |fc|>7V i=l i=0 

To prove the validity of the decomposition (13.21) it is enough to show that 



and 1 1 7m 1 1 < ^2- Taking into account (I3.48P and (I3.49p . the last means to give the estimate 



0<Ci< 



(3.51) 



This estimate is proved by Lemma 13.111 Therefore, the representation (13. 2p holds for any 
^ e M2, what completes the proof of the theorem. □ 



3.3 Auxiliary results 

In this subsection we prove several estimates which have been used in the proofs of Theo- 
rem 13.31 and Theorem 13. 4[ 

Lemma 3.5. Let us consider eigenvectors = ip{X^) and = 4'{X'^) and their represen- 
tations Ii2.18\) and ^2.19\) . Let us assume that ||a;J^||c" = ||z/mllc" = 1 in these representations. 
Then for any N E N there exists a constant C > such that 

ht\\<C, ^||^^||<C, m = l,...,£i, |fc|>iV. (3.52) 

In other words the two families of eigenvectors {yj,^ : m = 1, . . . ,ii,\k\ > N} with ||a^mllc" = 
1 and |p=^m • m = 1, . . . , ii,\k\ > with \\ym\\c" = 1 are bounded. 

Proof. Using (12.181) and the relation ||a;Jjj||c" = 1 we obtain: 

-1 

-1 

< 1 + e^RcA^p^if + < 1 + + 1^ < C, 

where r = maxr*^^) and r^'^^ are the radii of the circles Lt^{r^''^). The last inequality holds 

fceN 

since the real function -^'^ ^ decreases monotone from +cxd to 1 when the variable y runs 

y ^ 

from 00 to —0. 
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From fl2.19l) and since 



1, we have: 



A™ 



|2 Ill/mil^ + / 



-1 







d^ 



< 

< 
< 



..k ||2 / 1 



+ / 

-1 







.2ReAi 



p^P3(^)ll + 11^2(^)11) ds 



2ReA*L 



l+piFll^3WllL + P2 



L2 



+ 1) (P2WIIL + 1) + pipPsWIlL < c, 



L2(-i,0;C"X")- Here we used the fact that the real function ^^-^ 



where = 11^* 

increases monotone from to 1 when the variable y runs from 00 to —0. 

Remark 3.6. We notice that the norm of eigenvectors ip^ (assuming Hy^ll — ^) increases 
infinitely when A; — > 00. This could he seen on the example of eigenvectors of the operator 
A* (A;ie) = Am ^ 0): 



11^. 



k ||2 
mil 



k ||2 
mil 



+ 



-A^e ii2 

y-n 



-,2ReAS 



de = \\yLr\i + \x 



k |2' 



2ReA^ 



> (l + C|A^n ^ +00, k^oo. 
To formulate the next proposition we introduce the matrices 

/ ... 1 \ 



Rm 

I 



Rr, 



1 



. 

1 . 







1 
0/ 



m 



where / = In-m is the is the identity matrix of the dimension n — m. Obviously, Ri = I and 



R-} = R*= R for all m = 1, . . . 



Lemma 3.7. Assume that a i = {/ii, . . . , yU^^} consists of simple eigenvalues only. There 
exists N such that for any G Ai, \k\ > N and the corresponding matrix A(A,^) there 
exist matrices Pm,k, Qm,k of the form 



Pr 



( 1 -V2 

1 



m.k 



V 







-Pn \ 



1 J 



Qm,k 



( 1 

-^2 1 

V -qn 







(3.53) 



such that the product ji^Pm,kRm^{Xm)^mQm,k has the following form: 



A^ 

m 



Pm,kRm^{,Xrri)^'rnQm,k 



/ ... \ 







, det Sm,k 7^ 0, 



(3.54) 
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Moreover, for any e > there exists G Z such that for any > the components 
Pi = Pi{m, k), Qi = qi{m, k) of the matrices ^3. 531) may he estimated as follows: 



\Pi\ < hi <e, i = 2, . . . ,n. 
Proof. We begin with the analysis of the structure of the matrix 



(3.55) 



1 



Rm I A2{s) + jjrMs) I Rm ds 



Since the matrix y4_i is in Jordan form (12.161) . then the muhiphcation of A^i on Rm from 
the left and from the right changes the places of the one-dimensional Jordan blocks /ii and 



RmA^lRr, 

We introduce the notation 



S 



e^'Rm (Ao(s) + IAs(s) 1 Rm ds 



1 



A 



g g ^(n— l)x(n— 1) 



/ eii(A) ... ei„(A) \ 



\ £^ni(A) 



(3.56) 



.(A) J 



According to Proposition 14.61 elements of the matrix (13.561) tends to zero when |ImA| —>■ oo 
(and |ReA| < C < oo). Thus, |£^jj(A^)| — > when k oo. 

In the introduced notation the singular matrix ^i?mA(A^)i?m has the following form: 



—Rm'^{Xm)Rrn 



/ -1 +e~^™/^m + £ii(A^) ei2{\t) ... ein{\t 

^2l(A^) 



S, 



m,k 



where 



5 



m,k 



\ 



( 



\ 



£^nl(A, 



(3.57) 



/ £^22 (A 



'n-1 



+ e"-^'»5' + 



^2n 



\ ^n2(A^ 



Er, 



»(Am) / 



(3.58) 



and /„_i is the identity matrix of the dimension n — 1. Let us prove that 

det Sm,k 7^ 0. 



(3.59) 



Consider the identity —In-i + e ^"^S = —In-i + e ^"^S + (e — e ^"^)S, where Aj 
i(arg Hm + 27ik) is an eigenvalue of the operator A. Since e'^™ = /x^ "we have 



-I + e-^"^RrnA^iR„ 



^ 

-/„_i+e-^'"5 



and since the multiplicity of the eigenvalue £ o"! equals 1, we conclude that det(— + 



e-^-^) ^ 0. Since |Aj 



A 



such that for k : \k\ > N the estimates |e 



when k oo, then for any e > there exists > 
e-~H\S\\ < f and |e.,(A^)| < f hold. 
Thus, we have that det Sm,k = det(— J„_i +e~'^'"S' + Bm,k), where the absolute value of each 
component of -Bm,fc is less than e. Therefore, there exists A'" > such that 5*^,*; is invertible 
and we obtain the relation (I3.59p . 
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Since det Sm,k 7^ then from the relation (13 .5 71) we conclude that the first row of the 
matrix YrRm^{^m)^rn is a linear combination of all other rows and the first column is a 
linear combination of all other columns: 

^lj(-^m) = P2S2i + • • • + PnSni, i = 2,...,n 

^il(-^m) = Q2Sj2 + ■■■ + qnSjn, j = 2,...,n. (3.60) 

where Sij = Sij{m, k), 2 < i, j < n are the components of the matrix Sm,k- 

Let us consider the matrices Pm,k, Qm,k of the form (13.531) with the coefficients p2, . . . ,Pn 
and q2, . . . ,qn defined by (13. 3p . Direct computations gives us that YrPm kRm^{^m)RmQm k 
is of the form (13.541) . i.e.: 



/ ... \ 



in,k 



\ 







Let us estimate the coefficients p2, . . . ,Pn and q2, ■ ■ ■ , qn- The equations (13.31) may be 
rewritten in the form: 

Vl = {Sm^kfwi, V2 = Sm,k'W2, 

where Vi = (£i2(A^), . . . , ei„(A^))^, Wi = {p2,...,PnV, V2 = (£2i(A^), . . . , eni(A^))^, W2 = 
{q2, qnf- Since det Sm,k 7^ then 

^1 = {s;^]kfvi, W2 = s-\v2 

and since the values £^ij(Ajjj) are small then to show (I3.55P we have to prove that the estimate 
\\S-\\\ < C, C > holds for all k : \k\ > N. 

As we have shown above Smk = —In-i H — ^"S* + B^k, where elements of the matrices 
Bm,k tend to zero when k —>■ 00. Thus, there exists N E Tj such that for all : |fc| > 

^ dcf / \ ~ "'^ ^ 

the norm of the matrix Bmk = — { —In-i + — •S' ) Bmk is small enough, say ||-Bmfcll < k- 

Thus, there exists the inverse matrix of — B^^k for every \k\ > N, and these inverse 
matrices are bounded uniformly by k: 



Thus, we obtain the estimate 



J2{B^,ky\\ < Ci, \k\>N. 



I 9"^ I 

I m.k I 



{In- 



B. 



m,k ) 



i=0 



'n-1 



< c, 



what completes the proof of the lemma. 



□ 



Corollary 3. 8. Th e matrix function A„_fc(A) = jPm,kRm^{><)RmQm,k, where Pm,k, Qm,k 
are given by hS. 53|) . allows the following representation in a neighborhood U{\^) of the 
corresponding eigenvalue Aj^ G Ai, |A;| > A^.- 



^m,fc(A) 



/ (A-A^)r™''=(A) (A-A^)rr/(A) ... (X-XX'iX) \ 

: Sm,k{X) 

\ {x-xXi\x) ) 



(3.61) 
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where the functions r™'''^(A) are analytic in U{X^). Moreover, 



rn''(A^) 7^ 0, |r™''=(A: 



m,k/ \ k 



1, k 



oo. 



(3.62) 



Proof. Since A(A) is analytic, then all the components of the matrix function Am,A:(A) = 
jPm,kRm'^WRmQm,k are analytic in a neighborhood of the point A^. Moreover, since the 
matrix Am,fc(A^) has the form (13.541) . then we conclude that Am,fc(A) is of the form (13.611) . 

Let us prove the relation (I3.62p . If we assume that r^^' (Aj^) = 0, then (A — A^)r^' (A) = 
(A — A^)^f^''^(A), where r^i''{)\) is analytic. The last implies that the multiplicity of the 
root A = A^ of the equation det Am,fc(A) = is greater or equal to 2, i.e. det Am fc(A) = 
(A — A^)^r(A), where r(A) is an analytic function. Indeed, decomposing det Am,A;(A) by the 
elements of the first row, we see that all the term of this decomposition have the common 
multiplier (A — A^)^. Thus, we obtain that the multiplicity of A = A^ as the root of the 
equation det A(A) = is greater or equal to 2, what contradicts to the assumption that A^ 
is an eigenvalue of the multiplicity one of the operator A. 

Taking into account (I3.57P and the form of the transformations (I3.54p . we see that 

n n 

(A - Xt)rT{\X) = (-1 + e-Vm + en{X)) - $^P.^.i(A) - q,e^,{X). (3.63) 

1=2 j=2 

Differentiating (13.631) by A and substituting A = X'^, we obtain 



m,k 
^11 



(^m) = -e ^^/im + I eu{X) - ^PiEiiiX) -^^qjEijiX) 

i=2 



i=2 



The terms {Bij{X)y are of the form 



As 



sA2{s) + jMs) 



1 

A2 



ds. 



therefore, due to Proposition 14.61 and Lemma 13.71 we conclude that 

^0, k 



^11 (A) - ^PiSiiiX) -^qjEijiX) 

i=2 



oo. 



i=2 



A=Ai 



Since — e ^"^Hm — > —1 when k ^ oo, then we obtain the relation (13.620 and, in partic- 
ular, there exists a constant C > and an integer such that for \k\ > N we have 



< C < 



m,k/ \k 



11 



(A^) 



The last completes the proof of the proposition. 
Remark 3.9. The same arguments gives us that 



ra''(AjL)l-0, K-\Xl)\ 



m,k/ \ k 



OO. 



(3.64) 
□ 

(3.65) 



for all i = 2, . . . ,n and for all i 



n. 
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Proof. Indeed, let us consider r^''^(A) for j = 2, . . . ,n and use the fact that A_i is in a 
Jordan form: 

n 

(A - X'jrT/\X) = ei,{X) - J2p^e^AX) + Pj (-1 + e-V) +Pj-ic, 



2=2 



where /j, e a{A^i) and the constant c = if /i is geometrically simple or, otherwise, c = 1. 
Thus, we obtain 



j=2 



A=Ai 



□ 



and since pi = pi{m,k) when k oo due to Lemma I3.7[ then we conclude that 
|r^''^(Aj^)| when k ^ oo. 

Remark 3.10. Direct computations give us: 

/ 1 P2 ... Pn\ 
1 ... 



p-i 

m,k 



\ 



... 1 





/ 1 


. 






Qm,k 




1 . 


. 


(3.66) 




\ Qn 


. 







Lemma 3.11. Let o\ = {/ii, . . . , fie^^} consists of simple eigenvalues only. There exist con- 
stants < Ci < C2, N E Z such that for any A^ G Ai, \k\ > N the following estimate 
holds: 

^ ■{A'iXDxi^vt) 



< Ci < 



Ai 



<C2 



(3.67) 



where A'(A) = ^A(A); = x(A^), y'^ = y{X^) are defined by ^M), ^KW and = 

Wviw = 1- 

Proof. First, we prove the estimate fl3.67p for eigenvalues A^^ G Ai. Since G KerA(A^), 
then 

= ^p-,iPi,,A(At)Qi,feg->t = p-,iAi,,(At)Q->t, (3.68) 

where Pi,fc, Qi,fc, Ai,fc(A^) are defined by (13331) . (1331) . Thus, g~],x^ G KerAi,fc(A^) and, 
taking into account the form (13.541) of the matrix Ai_fc(A^), we conclude that Q^lx^ = 
(xi, 0, . . . , 0)"^, Xi 7^ 0. On the other hand, multiplying directly QTfc given by (I3.66P on 
x\ = ((x^)i, . . . , {xi)n)'^, we conclude that xi = {xi)i and (xj),; = —qi{x\)i, i = 2, . . . ,n. 

Due to the relation (13.551) . for any e > there exists G N such that for all k : \k\ > N 
we have: 

1 = \\4r = (1 + |g2p + . . . + |g„P) < \{x1),\' (1 + (r^ - 1)^') 

what implies that |(xj)i| —>■ 1 when k —>■ 00. Finally, we obtain 

Qr,M= ((^?)i>0,...,0)^, 0<C< |(xt)i| <1, \k\>N. 
Conjugating (I3.54p . we have the relation 
/ ... \ 



(3.69) 








\p,,,A{Xl)Q,A = 1=QI,A*{X1)PI, = AIM). (3.70) 



xt 



At 
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Let us use the fact that G KerA*(A^): 
A? 



(3.71) 



From the last we obtain that {Pi f.) Hi G KerA^;,(Aj), and, taking into account the left- 



hand side of f l3.70p . we conclude that (P*^ 



{yi,0,...,Of, yi ^ 0. Multiplying 



(Pi*fc) ^ ony'i = ((?/i)i, . . . , (yi) i)^ we obtain the relations yi = {y^)i and (?/J^)i = -Pi(?/i)i, 
i = 2, . . . ,n. Thus, due to fl3.55p . any e > and k : \k\ > N we have: 



fe||2 



(1 + Ip^r + • • ■ + IPnP) < (1 + - 1)^') 



and we conclude that \{yi] 



1 when k oo. Finally, 



{Pl^rVi = ((yi)i, 0, . . . , 0)^, 0<C< |(y^)i| < 1, \k\ > N. 
Differentiating f l3.6ip by A and putting A = Aj^, we obtain 



(3.72) 



/ rltiX'i) 



KM) = { jPi,kA{\)Qi,k 



A=A* 



^1,^ ( ^A'(At) - p|y^A(At) ) g,.. 



Using (13. 3p and the relation G KerA(A^), we obtain 



(3.73) 



^ (A'(At)x^, yf> = {Pi-^P,,, (3^A'(A^) - (x|FA(At)) gi,.gr>?, yi) 

(Pi,, (3:.A'(At) - (xfp A(At)) Qi,kQiixl {Pi^lTyl) (3.74) 



Finally, using the representation (13. 3p of the matrix /\\^{\\) and representations (I3.69p . 
(I3.72P of the vectors Qi\x'l, (Pf^ )*|/f, we conclude that 



A^ 



(A'(At)xt, y\) = r\^{\\){x\Uy\) 



(3.75) 



Moreover, taking into account the estimate (I3.62p of Corollary 13.81 and (I3.69p . (I3.72p . we 
obtain the estimate (I3.67p . what proves the lemma for the case of eigenvalues Af, i.e. for 
m = 1. 

Let us now prove the estimate for A^^ G Ai, m = 2, . . . , £i. In this case the idea of 

the proof remains the same but the arguing appears to be more cumbersome. In the proof 
we omit some detailed explanations that were given above for the case m = 1. 

Let us consider the product PmA(A^)Pm. Using the relation G KerA(Aj^J, we have 

= -^RrnPm,kPm,kRmA{^m)RmQm,kQm,kPmXrri — PmPm,kAm,k{^m)Qm,kPm'^m- (3.76) 
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Thus, Qmk-Rmx'^ G KerAm fc(A^) and from the exphcit form fl3.54p of Am,fc(A^) we conclude 
that Q'^^kRmX^ = (a^i, 0, . . . , 0)^, xi ^ 0. Multiplying Q'^^ on from the right, we 
changes places the first and the m-th column of therefore, we obtain: 



-qi[x. 



2, . . . ,n,i ^ m. 



Thus, taking into account (13.551) . for any e > there exists G N such that for all /c : \k\ > N 
we have: 



and, thus, \{x^)m\ 



1 when /c — > oo. Therefore, 

, T 



0<C<\{x'Jm\<l. \k\>N. 



The similar arguing gives us that 

{p-\rRmyt = {{yt)m, 0, • • • , 0)^, < C < \{ytU < l, \k\ > N. 
Differentiating (13.611) by A and putting A = Aj^, we obtain 



(3.77) 



(3.78) 



/ r-''{Xi) r^\\l) ... rr/(A^) \ 



ni,kf \ k 



21 V^m 



A' fA^ 



(A 



k \2 



^(^m) I RmQm,k- 



(3.79) 



Finally, using (13.771) . (I3.78p . (13. 3p and the relation G KerA(Ai), we obtain 



^ (A'(A^)xi, yt) = 

nPm,kPm,kRm i^m) ~ JJk-y^i^m)^ RmQm,kQm,kR"i^Tn^ V' 

^^m,ki.^ni)RmQ m,k-^ni-: {Pm,k) Rmyrr, 
^11 i^m)i-^m)m{ym.)m- 



To complete the proof of the lemma it remains to apply the estimates (I3.62p . (I3.77P and 



4 Boundedness of the resolvent on invariant subspaces 

In this section we prove the exponential stability of the restriction of the semigroup {e*'^}t>o 
onto (i.e. the semigroup {e^'^\M^}t>o)y where the invariant subspace is defined in 
Section 3 by (I3.33p . 

To show this we use the following well-known equivalent condition of exponential sta- 
bility (see e.g. [27, p.ll9] or [12, p.l39]): 

Let T{t) be a Co-semigroup on a Hilbert space H with a generator A. Then T{t) is exponen- 
tially stable if and only if the following conditions hold: 

1. {A : ReA > 0} C p{A); 

2. \\R{\A)\\ < M for all {A : ReA > 0} and for some constant M > 0. 
Thus, we reformulate our aim as follows. 
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Theorem 4.1 (on resolvent boundedness). Let cxi = {/ii, . . . , fiti} consists of simple eigen- 
values only. On the subspace defined by l\3. 3!^) the restriction of the resolvent R{X, A)\mo 
is uniformly bounded for A : ReA > 0, i.e. there exists C > such that 

||i?(A,^)x|| < xeM°. (4.80) 

Let us briefly describe the ideas of the proof. From the explicit form of the resol- 
vent (12.131) we conclude that the main difficulty is to prove the uniform boundedness of the 
term A~^{X)D{z,^, A) in neighborhoods of the eigenvalues of A located close to the imag- 
inary axis. Indeed, since det A(A^) = for A^ G Ai and ReA^ when k ^ oo then 
norm of A~^(A) grows infinitely when ReA — > and ImA — > oo simultaneously. However, 
the product A^^{X)D{z,C^, A) turns out to be bounded for {z,^{-)) G Mj. 

Lemma 4.2. The vector- function A^^{X)D{z,C,, X) '■ M2 x C" ^ C" is uniformly bounded 
in neighborhoods Us{X^) of eigenvalues A^ G Ai for some fixed 6 > 0, i.e.: 

1. for any A; : |A;| > and m = there exists a constant Cm,k such that the 
estimate \\A-^{X)D{z,i,X)\\ < CmA\i.z^ii.-))hh holds for all X G f/5(A^) and{z,^{-)) G M^. 

2. there exists a constant C > such that Cm,k < C for all m = 1, . . . k : \k\ > N. 

The proof of this lemma is technically difficult. It essentially uses the following relation. 

Lemma 4.3. For any vector g = {z,^{-))'^ G and for any eigenvalue A^ G Ai the 
following relation holds: 

D{z,^,X'jelmA{Xt). (4.81) 
The complete proofs of the mentioned propositions are given in the next subsection. 



Remark 4.4. Theorem "JH holds also for the subspace defined by under the 



assumption A2 = 0. The proof remains the same. 

4.1 The proof of Theorem 14JJ 

We begin with several auxiliary propositions. 

Proposition 4.5. // the vector y G 1mA, A G C"^", then for any two matrices P, Q such 
that det Q the relation Py G lm{PAQ) holds. 

Proof. The relation y G 1mA means that there exists a vector x such that Ax = y. Since Q 
is non-singular then there exists a vector xi such that x = Qxi. Therefore, AQxi = y and, 
multiplying on P from the left we obtain PAQxi = Py. □ 

Proposition 4.6. Let Lq G C be a bounded closed set and f{s) G L2[— 1,0]. Denote by 


(^k{X) = f e^'^^'^^e^'^ f{s) ds, X G Lq, k G Z. Then |afc(A)| — > when /c — > 00 uniformly on the 

0- 



-1 

set Ln. 



Proof. Integrals afc(A) can be considered as Fourier coefficients of the function e'*'^/(s), thus, 
they converge to zero when k ^ 00. It remains to prove that they converge uniformly on 
the set Lq. The last means that for any e > there exists n G N such that for any |A;| > n 
and for any A G Lq we have |afc(A)| < e. 

Let us suppose the contrary: 3e > such that Vn G N, 3|A;| > n, 3A G Lq: |afc(A)| > e. 
Thus, there exists a sequence ki < k2 < ■ ■ ■ and a sequence {Xk-}°^^ such that \ak-{Xk-)\ > e. 
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Since Lq is a bounded set then there exists a converging subsequence of {A^-}^^ which 
we denote by {Xj}j^j, where J C is a strictly increasing sequence. Moreover, since Lq is 
also closed, then the limit of belongs to Lq: \j — > Aq G Lq. Let us show that the 

sequence {ajt(Ao)} does not converge to zero. 

Indeed, choosing big enough n G N, such that for any j > n, j E J and any s G [—1,0]: 



e ^■''1 < I II /(s) II, we obtain 



|aj(Ao) - aj{Xj)\ 



< 



e'^Vis) ds<-. 



Since |aj(Aj)| > e and assuming that |aj(Ao)| < |aj(Aj)|, we obtain 



2 > Wji^o) - aj{^j)\ > \aj{^j)\ - \aj{^o)\ >e - |%(Ao)|, 

and, thus, |aj(Ao)| > | for any j E J, j > n. 

Thus, {flfcl-^o)} does not converge to zero and we have obtained a contradiction with 
the fact that they are the coefficients of the Fourier series of the function e^°*/(s). The last 
completes the proof of the proposition. □ 



oo and — oo < a < ReAz, < b < oo 



Corollary 4.7. If the sequence {Xk} is such that ImA^ 
then for any f(s) G -^^2(0, 1; C"^™) we have: J^^ e'^^^f(s) ds when k ^ 00. 

Lemma 4.8. The following estimates hold: 

1. There exists a constant C > such that ||iA(A)|| < C for all A G {A : ReA > 
0}\f/(0), and ||A(A)|| < C for all A G f/(0). 

2. There exists a constant C > such that || ^, A) || < C*! (z, ,^(-)) Ha/s for all 
A G {A : ReA > 0}\f/(0), and \\D{z,i,\)\\ < C||(^,e(-))llM, for all A G U{Q), {z,i{-)) G M^. 

Proof. From the explicit form fl2.14p of A (A) we have an estimate 



A 



A(A) 



< 1 



l^-i|| + 



j e^M2(s)ds y e^M3(s)ds 



for A G {A : ReA > 0}\f/(0), and, thus, it remains to prove that there exists a constant 



Ci > such that 



/_V^^.(s)ds 



< Ci, i = 2,3. Indeed, if we suppose the contrary, then 



there exists an unbounded sequence {Xj}f=i such that \\j\e^^'Ai{s)<ls 



00 when j 00. 



On the other hand, it is easy to see that for any A; > 0: J_-^ e'^^s'^ds — > when |A| 00 
and A G {A : ReA > 0}. Since the set of polynomials is everywhere dense in L2(— 1,0), 



then 



/_V^A(3)ds 



when |A| — 00, A G {A : ReA > 0} and we have come to a 

contradiction. 

The estimate ||A(A)|| < C for all A G f/(0) follows easily from the explicit form fl2.14p . 
The estimates for ^, A) may be checked directly in the same manner, taking into 
account that e~^ e~^''s''ds when |A| — > 00 and A G {A : ReA > 0}, > 0. □ 



Now we pass over to the proofs of the main propositions mentioned in the beginning of 
the section. 

Proof of Lemma \4.2\ Let us introduce the following notation: 



fiXf^A~\X)Diz,C,X) 



-A(A) 



D{z,^,X)), {z,a-))eMl 



(4.82) 
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We analyze the behavior of the vector-function /(A) near the imaginary axis. For the points 
e Ai, which are the eigenvalues of the operator A, the inverse to the matrix A(Aj^) does 
not exists. These eigenvalues approach to the imaginary axis when A; — > oo. Our first aim is 
to prove that /(A) is bounded in each neighborhood f/(A^) of A^ G Ai, i.e. that the limit 
lim A-\\)D{z,^,\) exists for all (^,{(-)) e Ml 

Since A(A) and D{z,C,,\) are analytic and since, by the construction, all eigenvalues 
A^ G Ai are simple, then we have that if A^ is a pole of /(A) then it is a simple pole. In 
other words, in every neighborhood U{\^) the vector-function /(A) may be represented as 
follows: 

fW = ^^/-l + - ^'rnYU (4.83) 

'■ i=0 



Thus, our aim is to prove that for each A^ the coefficient / 
to zero in the representation fl4.83p . i.e. that /(A) is analytic. To prove this, we construct a 



lim (A — A^)/(A) is equal 



representation of the matrix (iA(A)) ^ which separates the singularity of this matrix. 

According to Lemma 13. 5^ for each A^ G Ai there exist matrices Pm,k, Qm,k such that 
the value of the matrix-function Am,fc(A) = jPm,kRmAWRmQm,k at the point A = A.^ has 
the form fl3.54p . i.e. 



/ 



Am,fc(A, 



'Pra,kRmA{^^rn}R"i^Qi^yh 



\ 



m,k 



\ 



det Sm^k ^ 0. 



We rewrite the representation fl4.82p of the function /(A) in a neighborhood U{\^ as follows: 
/(A) = (ii?„P^iP„,fci?„A(A)i?„Q„,fcQ-|,i?„)"'(iZ}(z,e,A)) 

= RmQm,k {jPm,kRmAWRmQm,k) Pm,kRm {jD {z , ^ , \)^ (4.84) 
= RmQm,k [^m,kW) Pm,kRm {jD , , ^)) ■ 



Let US consider the Taylor expansion of the analytic matrix-function Am,A:(A) in f/(Af„): 

oo _ 

A„,.(A) = A^A^t) + (A - A^)A'„,,(A^) + E 7t(^ " A^)^^2^(A^)- (4-85) 



i=2 



Due to Corollary (13. 8p . Am,kW allows the representation (I3.6ip in some ?7(A^), i.e 

/ (A - A^)r-''=(A) (A - A,i)rr/(A) ... (A - A^)rr/(A) \ 
(A - A^)r-''(A) 

: 'S'm,fc(A) 
V (A-A^)C^(A) 



Am,fc(A) 



where ^■^''^(A) are analytic functions, and we note that Sm,k{^m) 
defined by (13.540 . Differentiating the last relation by A at A = A^, we obtain: 



Sm..k, where S^ h is 



/ ---'(A^) rT,\Xi) ... rr/(A^) \ 



^m,k(^r. 



m,k/ \ k 



21 V^m 



\ ^nl (Am) 



^'m,ki^rn) 



m,fc ' m,ki 
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■^0 def 

m,k 











V 







m,k 



/ ^0 

^21 (Am) 



... \ 



'-^m,k\ ml 



V ^nl (Am) 



We introduce the matrix-function Fm,jt(A) =^Am,fe(A^) + (A — A^)r^ j^, which has the following 
structure: 

f rT{\Xi){X - Xi) rT,\\t)i\-\i) ... rr/(A^)(A-A^) \ 


; Sm,k 





-^m,fc(A) 



(4.86) 



/ 



The matrix Fm,A;(A) is non-singular in a neighborhood U (Aj^)\{Aj^}. Indeed, due to Lemma 13751 
and Corollary 13.81 we have that det S'm,fc 7^ 0, r^' (A^) ^ 0, and, thus 

detFm,fc(A) = T^\\l){\ - A^)det^m,fc 7^ 0, A G f/(A^)\{A,^}. 
Therefore, there exists the inverse matrix i^^\,(A), which is of the following form: 



^m,fc(-^) 



1 



rpTn,k 

r-n''(A^)(A-A^) ^21 

F^/ 



-:m,k 
'2n 



T-im,k \ 
rpm,k 



'<in,k 



(4.87) 



where 



T-im,k 



1 ' ' ^ k k 

r7]-^(A^.) dot ^J^^'^y^^^lj iK)[Sm,k{K)]ij, 



i = 2, 
2, 



(4. 



and by [Sm,kW]ij we denote the complementary minor of the element s^'^(A), i,j = 2,.. 
of the matrix ^^.^(A). Since the matrix- functions S'm,fc(A) are analytic and since 5'm,fc(Aj^ 



, n 



S when k ^ oo, then ||S'm_fc(A)||, || [5'm,fc(A)]jj|| and \s^' (A)| are uniformly bounded for all k 
and A e Us{X^). Thus, we conclude that l-Fj^ '^l — ^ ^ ^! J = 2, . . . , ra. Moreover, 



smce r 



m,k/^k 



(AJ 



due to Remark 13.91 then F^'^ ^ 0, z = 2, . . . , n when k — + oo. 



Let us rewrite the representation (14.851) as follows: 

/^m,kW = i^™..(A) + (A - A^^)r^ + E - ^™)^^!:y^™ 



i=2 

Fm,.(A) ( / + (A - A^)F„i(A)r^ + E - ^typ^^^^kW^ZM 

i=2 



and introduce the notation 

Tm,fc(A) ='^(A - A^)F~j,(A)r;^ + ^ -(A - A^)*F~\.(A)A^yA, 



(4.89) 



i=2 



Let us prove that for any e > there exist 6 > and G N such that for any k : |A;| > 
the following estimate holds: 



|Tm,.(A)||<e, XeUsiX'J. 



(4.90) 
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From (14.891) we have that jt(A) = F^\{X)Am,k{^)—I, and we are proving the estimate 
||i^^yA)A^,fc(A) -I\\<e, Xe Us(Xt)- Using the representations (EH]), fOTl) and (Km . 
we estimate the elements {7ij''^(A)}"j=i of the matrix Tm,kW- For the sake of convenience, 
we divide these elements onto several groups and we begin with the element 7ii'^(A): 

m,k / \ \ n 
"^11 \^m) i=2 

Due to Corollary 13.81 and since r^i^[X) is analytic, there exist 5 > and G N such that 



m.K f \ \ 



„ < ^ for all k : \k\ > N and A G Us{X^. Besides, since F^'^ —>■ 0, i = 2, . . . ,n 

^11 \'^m) 

and r^'''{X^) — when A; — > oo, we obtain the estimate 7n'^(A) < ^ for all /c : \k\ > N, 
X G t/5(A^). 

Let us consider other diagonal elements of the matrix T„ifc(A): 

n n 

iTi'W = E Fi;''sTi\X) - 1 = 5:(-l)^-^^[5^,.(A)].,(.-'^(A) - sTi\Xt)), J=2,...,n. 

i=2 1=2 

There exists S > such that |7j|J''^(A)| < ^ for all k : > A^, A G ^/^(A,^). Further, we 
consider the elements of the first row: 

m,k / \ \ n 
m,k / w Ij \ ' \ ^ T-,m,k m,k / \\ ■ o 

Since F™''' — ;> 0, i = 2, ...,n and ^"^''^(A^) when A; ^ oo, we obtain the estimate 

7Tj''W\ < t for all k: \k\>N,Xe Us{Xl). 
Finally we consider all other elements: 

n n 
r=2 r=2 

They may be estimated as |7j'J''^(A)| < ^ for all k : |A;| > A^, A G Us{X^) choosing small 
enough 5 > 0. 

Finally, we obtain the estimate (14.901) and, therefore, there exist 5 > 0, G N such 
that the matrix / + T^ k^X) has an inverse for any A G ^/^(A^), k : \k\ > N: 

(/ + T^,fc(A))-i = / + (A - A^)r„,fc(A), (4.91) 

where rm,fc(A) is analytic in a neighborhood Us{X^). 
Finally, from KMf . fICTD and fOTj) we obtain: 

/(A) = RrnQm,k^-^^k{^)Pm,kRm{\D{z,i,X)) 

= RmQm,k {F^,kW{I + T^,fc(A)))-' RmPn.,k {{D{z, ^, A)) (4.92) 

= RmQm,k (/ + (A - Xi)Tm,kW) F-\{X)P^,kRm {\D{z, ^, A)) . 

Let us use the fact that {z,C,{-)) G Mg. In Lemma [4.31 the important relation D{z,^, X'^) G 
ImA(A^), {z, (,{■)) G M2 is stated. Thus, due to Proposition 14.51 we conclude that 

^Pm,kRmD{z, A^) G ImA„,fc(Aj;,). (4.93) 
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Moreover, since the matrix A,„ fc(A^) is of the form fl3.54l) . we conclude that the first 
component of the vector -rrPm kRmD{z, ^, X^) equals to zero: 

^Pm,kRmD{z, e, A^) = dm,k = (0, C2, . . . , C^f , (4.94) 

and since the vector-function jPjn,kRmD{z,C,, X) is analytic in a neighborhood U^X'^), we 
conclude that ^ 

jPm,kRmD{z, A) = dm,k + rfm,fc(A), rf„,fc(A^) = 0. (4.95) 

Finally, we note that F^\{X)dm,k is a constant vector and the vector-function F~\{X)dm,k{X) 
is bounded in U{X'^). Taking into account fl4.92p . (14.871) and (I4.95p . we obtain that 

hm (A - A^)/(A) = R^Q^,, lim (A - A^)F„-i(A)P^,fci?„ (iD(^, A)) 

= i?^Q„,fc limJA-A^)F-i(A)(d;,,fc + d„,fc(A)) (4.96) 
= 0. 

Thus, we have proved that /(A) = A~^{X)D{z,^, X) is an analytic vector-function in 
Us{X^), k : \k\ > N, m = 1, . . . ,ii what gives the estimate \\A~^{X)D{z,^, X)\\ < Cmk, 
iz,a-))^Ml 

It remains to prove that /(A) is uniformly bounded in the neighborhoods Us{X'^) for all 
k : \k\ > N, m = 1, . . . In other words, our next aim is to prove that the set of vectors 

/o = f^'' = f{Xt) = {A'\X)D{z,^,X)),=,.^ 

is bounded. Taking into account the representation (14.921) . we obtain: 



fr = {Rn.Qn.,kF;,l{X)Pm,kRm{\D{z,^,X)))^_^, 
= RmQm,kFm.kWid'm,k + dm,k{X)) \=\k^ 



1 l^m; j=2 i=2 



m,k 



(4.97) 



where c?]^ = (c^mfcl-^m))! the first component of the derivative of dm,k{X) at the point A^. 

As we have mentioned above, there exists a constant Ci > such that < Ci, 

for all A; G N, m = 1, . . . , £i. The estimates ||-Pm,fc|| < C*! and ||Qm,A;|| < C*! follow from the 
estimate (13.551) of Lemma 13. 7[ The estimates | Cj | < Ci and \d\\ < Ci follow immediately 
from Lemma 14. 8[ From the relation (I3.62p of the Corollary 13.81 it follows that there exists a 
constant C2 > such that < C2 < |'"ii''^(A^)| for all A; G N, m = 1, . . . , di and, thus, 

1 1 
— < — . 



11 V^m 



Finally, we conclude that ||/™''^|| < C for all m = 1, . . . , £1, /c : \k\ > N, what completes 
the proof of the lemma. □ 

Proof of Lemma \4.3\ Since g G M° then c/±Vm = V'(A^) for all A^ G Ai. Therefore, the 
proposition follows from Lemma [2.91 □ 
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Proof of Theorem \4-l\ Let 5 > be such that Lemma [4.21 holds. We divide the closed right 
half-plane onto the following two sets: 



K,{6) = {A : ReA > 0, A G Us{\t), Xt e 
K2{6) = {A : ReA > 0}\Ki. 



(4.98) 



First, let us estimate ||i?(A,^)x|| for any A G Ki{S) and x G M^. Due to Lemma [4.21 
we have: \\A-\X)D{z,^, X)\\ < Ci\\x\\, x = iz,^{-)) G M^. Due to Corollary we have 
the estimate || J^ie~'^^^is) ds\\ < C2||x||. Thus, for any x = {z,^{-)) G Mg, A G Ki{5) we 
obtain: 



mx,A)x\ 



e-M_i / e-^'^{s) ds + (/ - e-^A_i)A-\X)D{z, ^, A) 



+ 



as) ds + e''A-\X)D{z,^,X) 



< e^||A_i||C2||x|| + (1 + eiA_i||)Ci||x| 



L2 



(4.99) 



< 



+ / / e^^'-'^^s) ds + e^'A-\X)D{z, e, A) 
V-i 

e^\\A_i\\C2 + (1 + eiA_i||)Ci + (e^Ca + Cf)^ 



d^ 



Let us consider A G K2{6). There exists e > such that |^det A(A)| > e for any 
A G K2{5)\U{Gi). Indeed, if we suppose the contrary then there exists a sequence {Xi}^^ 

y-det A(Aj) — s> 0, i — » oo. If the sequence {Xi}^^ is bounded, then it possesses 

0. However, the closure of 



such that 



idet A(A) 

A 



a converging subsequence: Xi. X and, thus 

the set K2{5)\U{Q) does not contain zeroes of the function det A(A) and we have obtained 
a contradiction. 

If the sequence {Aj}^^ is unbounded, i.e. |Ai| oo when z ^ oo, then we have 
j^^e^''^ A2{s) ds ^ and ^ j'^^e^^'^ A'i^s) ds 0, i ^ oo. Moreover, since | det(— / + 
e~^A_i)\ = I Yli^ + e~'^/im)| ^ n |1 + ^Vml, we conclude that 



A 



det A (A, 



det(-/ + e-^^A^i + J e^^'A2{s) ds + — y e^'*A3(s) ds 74 0, i 00. 



Thus, we have obtained a contradiction again. 

Taking into account the estimates ||;|A(A)|| < Ci and ||^D(z,,^, A)|| < C2||x|| from 
Lemma SSI we conclude that \\A-\X)D{z,^, X)\\ < CsWxW for all A G K2i5)\U{0). It is 
easy to see that ||e"^ e~^''^(s) ds|| < C4||x|| for all A G K2{6)\U{0). Finally, similarly 
to (14.991) we obtain the following estimate 



\\R{X,A)x\\ <C\\x\\, XeK2. 
The last completes the proof of the theorem. 



□ 



5 Stability analysis 

Basing on the results from Section [3] and Section HI we prove the main result on stability 
which does not assume the condition det A_i 7^ 0. 
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Theorem 5.1 (on stability). If (y{A) C {A : ReA < 0} and Ox = o-{A_i) n {n : \fi\ = 1} 
consists only of simple eigenvalues, then system ( [i. ?| j is strongly asymptotically stable. 

Proof. Let us show that ||e*-^x|| when t +oo for any x G M2. Due to Theorem 13.31 
each X e M2 aUows the foUowing representation: 

X = Xq + Xi, Xo e M2, Xi E M2 , 

where and are defined by fl3.3ip - fl3.33p Moreover, the basis of consists of the 
following eigenvectors: 

{^l : {A - = 0, A^ G Ai = Ai(iV)}. (5.100) 

Thus, for any xi G M\ we have the representation 

^1 ^1 ^1 

m=l |fc|>Af m=l |fc|>iV m=l |A;|>Af 

Let us consider a norm || ■ ||i in which the Riesz basis (15.1001) is orthogonal, then we have 
the following estimate: 



m=l |fc|>Ar 

Since the series ^ ^ cJ^V^m converges and since Hv^mlli < C for all k and m = 

m=l \k\>N 

h 2 
l,...,£i, then for any £ > there exists Ni > N such that ^ ^ l|CmV^mlli — 

m=l |fc|>Ari 

Moreover, since the set {(m, fc) : m = 1, . . . , £1, < < ^"1} is finite and since ReA^ < 0, 
then there exists to > such that for any t > to we have: ^ ^ 6^^'''^'"*||c^(y9^||^ < ^. 

m=l Ar<|A:|<Ari 

Thus, we obtain 

E E < E E -'^^^'14^^11? + E E W'tvtwi < j- 

m=l |A:|>Ar ■m=l Ar<|fc|<Ari m=l |fc|>Afi 

(5.102) 

Due to Theorem 14. II the semigroup e^-^ljv^o is exponentially stable, i.e. by definition there 
exist some positive constants M, lj such that ||e*'^|^.fO || < Me""^*. Thus, for any xo G M2 
there exists to > such that for any t > to we have an estimate 

||e*^a;o||i < Me-'^*||xo||i < |. (5.103) 

Finally, from the estimates fl5.10ip , (15.1021) and fl5.103l) we conclude that for any x G M2 
and for any e > there exists to > such that for any t > to the following estimate holds: 

lle^-^xlli < ||e*-^xo||i + ||e*-^xi||i < e. (5.104) 

Therefore, lim ||e*'^x||i = 0, what implies that the system (11. 6p is strongly asymptoti- 

t— >+oo 

cally stable. □ 
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5.1 An example of dilemma: stable and unstable situations 

In this subsection we give an explicit example illustrating the item (iii) of Theorem II. 1[ 
Namely, we construct two systems having the same spectrum and satisfying the following 
conditions: cr{A) C {A : ReA < 0} and there are no Jordan blocks, corresponding to 
eigenvalues from (Xi = cr(yl_i) fl {/i : = 1}, but there exists an eigenvalue ^ E cxi whose 
eigenspace is at least two dimensional. Moreover, one of the constructed systems appears to 
be stable while the other is unstable. 
We consider the system of the form 

iW=(~J zeC',t>0, (5.105) 

where 6 is a real positive number and for the value of s we essentially distinguish two cases: 
s = and s ^ 0. 

Remark 5.2. The systems of the form z = A_iz{t — 1) + AQz{t) is a special case of the 
systems ( [i. T] ). 

Proof. Indeed, we choose ^2(6*) = {9 + 1)Aq and ^3(6') = Aq. Thus, 

/O pO pO nO 

A2ie)z{t + 9) d9 + J A3{9)z{t + 9) d9 = J (9 + l)Aoz{t + 9) d9 + J Aoz{t + 9) d9 

= Aoj {{9 + l)z{t + 9)y d9 = Aoz{t). 

□ 

The eigenvalues of the operator A are the roots of the equation det A^(A) = 0, which, 
in our particular case, has the form: 

det(-A/ + Xe-^A_i + Aq) = det ^ " ^ - ^ ^ _ 6 ) = 0- 

Thus, all the eigenvalues of the operator A satisfy the equation 

Ae^ + A + be^ = (5.106) 

and the multiplicity of each eigenvalue equals two. To prove that (t{A) C {A : ReA < 0}, we 
use the results on transcendental equations obtained by L. Pontryagin [17]. 

Let us consider the equation H{z) = 0, where H{z) = h{z,e^) is a polynomial with 
respect to z and e^. 

Definition 5.3. We say that the function H{z) = ^ amnz"^^^^ possesses the principal term 

m,n 

ttrsZ^Q^^ , if for all other terms am„z™e"^ we have that r > m and s > n. 

We denote by F{y) : M ^ M and G{y) : M — > M, correspondingly the real and the 
imaginary parts of the function H{iy), i.e. H{iy) = F{y) + iG{y), y eM.. 

Definition 5.4. We say that the zeroes of two real-valued functions of a real variable alter- 
nate if and only if 

a) each of these functions has no multiple roots; 

b) between any two zeroes of one of these functions there exists at least one zero of the other; 

c) the functions are never simultaneously zero. 
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Theorem 5.5 ([T71 Pontryagin, 1942]). Let H{z) = h{z,e^) be a polynomial with a principal 
term. 

1. If all the zeroes of the function H{z) belong to the open left half-plane: ReAfe < 0, then 
the zeroes of F{y) and G{y) are real, alternating and for ally (zM. the following inequality 
holds 

G'{y)F{y)-G{y)F'{y)>Q. (5.107) 

2. Any of the conditions below is sufficient for all the zeroes of the function H{z) to lie in 
the open left half-plane: 

a) all the zeroes of the functions F{y) and G{y) are real, alternating and the inequality 
l^5.101\) is satisfied for at least one value of y; 

b) all the zeroes of the function F{y) are real and for each zero y = yo the inequality 
li5.1U7\ ) is satisfied, i.e. G{yo)F'{yo) < 0; 

c) all the zeroes of the function G{y) are real and for each zero y = yo the inequality 
^5.1U7\ ) IS satisfied, i.e. G'{yo)F{yo) > 0. 

The following theorem gives a criterion for all zeroes of a function to be real. 

Theorem 5.6 ( JTli Pontryagin, 1942]). Let F{z) = f{z, cos z, sin z) be a polynomial with a 
principal term z^(j)i^{cosz,smz), where {cos z,sm z) is homogeneous with respect to cos 2 
and sin z polynomial. 

The function F{z), 2; G C possesses only real zeroes if and only if for all big enough A; G Z 
the function F{x), x G M possesses exactly Aks + r real roots on the interval —27ik -\- e < 
X < 2'Kk + e for some e > 0. 

We use the results mentioned above to analyze the location of the roots of the equa- 
tion (ElonD. 

Proposition 5.7. For any b > and any s G C the spectrum of the corresponding operator 
A belongs to the open left half-plane. 

Proof. We represent H{iy) = iye^^ + iy + 6e*^ as follows: 

H{iy) = iy{cos y -\- i sin y) -\- iy -\- b {cosy -\- ism y) 
= {b cosy — y sin y) -\- i{y cosy -\- y -\- sin y) 
= F{y)+iG{y). 

Since the roots of the equations ?/ = and cosy = are not roots of the equation F{y) = 0, 
we rewrite F{y) = as tgy = ^. Several zeroes of the equation F{y) = may be seen on 
the figure below: 
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I I III I \i\ I I 
-2.5 





x=b/y 



Figure 3. 

We see that exactly 5 zeroes of the equation F{y) = belong to the interval —27i + e < 
y < 2tt + e if we choose e > " big" enough such that the zero on the right from 27r belongs 
to the interval. Adding to this interval from the left and from the right 27r we add each time 
4 other zeroes. Thus, on each interval —2TTk + e < y < 2TTk + e we have exactly 4k +1 zeroes 
of the equation F{y) = and, therefore, the conditions of the Theorem 15.61 are satisfied. 
Thus, all the zeroes of the equation F{y) = are real. 

Let us now prove that for each root yo of the equation F{y) = the inequality 
G{yo)F'{yo) < holds. We use the notation cosy = C, siny = S: 



-GF' 



-hS-S- yoC){yoC + yo + hS) = {hS + S + yoC){yoC + + hS). 



We rewrite tgyo = — as C = — and substitute C at the last relation: 



-GF' 



bS + ^S + s] (^S + bS + yo 



Since > and 6 + ^ + 1 > 0, then it remains to prove, that the third multiplier is greater 

b 



then zero. From tgyo ~ ^ ^^^'^ have b 



1 + 6 + ^ 
S 



VqC , yoS yo 
S 



S 



c 



2/0- 



^ and substituting we obtain: 
l + G 



c 

C^ + S^ + C 



sc 



Va- 



se 



1 + C)bC' > 0. 



We note also that the fact that yo'S'C > can be easily seen from the picture above, 
and use the sufficient condition b) of Theorem 15.61 

Applying Theorem 15.61 we complete the proof of the fact that all the eigenvalues of the 
operator A belongs to the open left half-plane. □ 

Proposition 5.8. If s = then the operator A possesses eigenvectors only, i.e. it possesses 
no root vectors; if s then to each eigenvalue X G cr{A) there corresponds a pair of 
eigenvector and root vector of the operator A. 
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Proof. If s = then A_4(Afc) is the zero-matrix for each eigenvalue A^, and, therefore, the the 
space of solutions of the equation A_A^{\k)z = is two-dimensional. We choose the following 
basis of this space: zl — (1, 0)^, zl — (0, 1)^. 

Let us consider the equation for eigenvectors: 

The solution of the second equation is given by z{9) = e^''^z{0), thus, z{—l) = e~^''z{0). 
Taking into account the domain of the operator ^: y — z{0)—A_iz{—l) = {I—e~^''A_i)z{0), 
we obtain from the first equation: {Aq — X^I + Xke~^'' A_i)z{0) — 0, or, in an equivalent form: 
A^(Afe)^(0) = 0. 

As we have noted above the last equation has two-dimensional solution: z{0) = zl = 
(1,0)^ and z{0) = zl = (0, 1)"^, and therefore, there is a two-dimensional eigenspace of the 
operator A corresponding to the eigenvalue A^. Two eigenvectors of this subspace can be 

chosen in the following form: fl = ( if \ \ where vl = \ ^ ^ \ ^li^) ~ ^ 



= ( J-) ) ' ^ ( 1 + e-^^ ) ' "^^^^ ^ ( e° ^ ) • 

Thus, to any eigenvalue A^ of the operator A there corresponds the two-dimensional 
eigenspace. 

e"^*= 



If s = 1 (or s 7^ 0) we have that A_4(Afe) ~ q y space of 

solutions of the equation Ayi{Xk)z = is one-dimensional: zl = (1,0)^ (since, obviously, 
g-Afc _^ Thus, the equation A_A^{Xk)z = zl has also one-dimensional solution which we 
denote by zl = (0, 1)^. 

Arguing as above we show that the operator A possesses one eigenvector corresponding 

to the eigenvalue A^: fl = ^f^^^ ^ , where = ^ ^ + J ^ , zl{0) = ) ' 

Let us show that the operator A possesses one root vector. Each^root vector / of 
the operator A satisfies the following equation: {A — Xkl^f — 0, i.e. f = {A — Xkl)f G 

Ker(^ — Xkl). Prom the last we conclude that f — fl = i ^^x^e g^r 1 the root 
vector satisfies the relation 

Aoz{-l)-Xky \ _f (1 + e-^SOf 
m-hzie) )-[ (e^^^0)^■ 

The solution of the second equation is given by 

re 

z{e) = e^^^z{Q) + / e^^^^-^^e^^^yl dr = e^'=^^(0) + Oe^^^yl, 
Jo 

what implies z{—l) = e^^''z{0) — e~^''yl. Taking into account the domain of the operator 
A, we write down the first equation in the following form: 

^o^(O) - Xk {z{0) - A_,{e-^'^z{0) - e-^'^yl)) = (1 + e-^'^)yl 

which is equivalent to 

1 + e-^" - Xke-^>' 



A^(Afe)z(0) = (1 + 6-^*=)^^ + XkA^ie-^^yl 
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Thus, if 1 + e — Xf^e ^ then z{0) is the root vector of the matrix A_4(Afe) and 

yl 



therefore, z{0) = zl = (0, 1)^ and the operator A possesses the root vector f^=i 



where yl=y^^ j and zKO) = ^^^^ 

It remains to show that 1 + e"'^* — A^e"'^* ^ 0. We suppose the contrary: Xk^^^'' = 
1 + e~'^'= and, multiplying the last expression onto e'^*, we obtain = e'^'' + 1. Since A^ is an 
eigenvalue of A, then we obtain e^'=(e^'' + 1) + e^*" + 1 + e^'= = or e^^'' + 3e^'' + 1 = 0, and 
we conclude that e^'= = For the root e'^'= = "'^"^ we have that ReA^ > 0, and for the 

root e^*^ = we have: A^ = + 1 = ^^^^^ > 0. We have obtained the contradiction 

what completes the analysis of eigenvectors and root vectors of the operator A. □ 

Proposition 5.9. If s = then the system /[5.105\) is strongly asymptotically stable; if s 
then the system is unstable. 

Proof. First we consider the system flS.lOSp when s = and prove that in this case the 
system is stable. Let us evaluate the norm of eigenvectors. 

/O rO 
eA.V^e de=\l+ e-^*p + J e^i^^^^-^ d^ 

zReAfc 

Since ReAt. when k oo, then lim (l — 6"^^'^'^'=) = 1. Therefore, < Ci < 

(l — e"^^'''^'') < 6*2. Taking into account that |e~'^'=| < C3, we obtain the following 
estimates: 

Ci < Wflf < (1 + Csf + C2 = C4. (5.108) 

As it has been shown in ^Ij, the subspaces V^^'^^ = Lin{/jJ,/|} (and the finite-dimensional 
subspace W^) form a Riesz basis of the space M2. Since we have proved the estimate 
fl5.108p then the eigenvectors {/^, /|} (together with vectors from Wn) form a basis of M2. 
Therefore, we have a Riesz basis of eigenvectors. Further we just repeat the proof given in 
pn Theorem 23]. 

We consider a norm || ■ ||i in which the eigenvectors {/fc,/|}fcez orthogonal. Let a 
vector X belongs to a closed span of the subspaces V^^\ then x = {(^kfl + hfl) ^'^d we 



have: 



Therefore, 



e^'x = Y,^'''\ocufl+Ml). 



fcez feez 

and, thus, the family e'^* is uniformly bounded in the subspace generated by the subspaces 
V^''\ From the last we conclude that the system is strongly asymptotically stable. 
Let us consider the system (15.1051) when s = 1. 

Let an operator A has a sequence of eigenvalues {Afcj^i such that ReA^ < and 
ReAfc when k ^ 00 and to each Xk there corresponds one eigenvector f ^ and at least 
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one root vector . We show that the equation x = Ax is unstable. Let us suppose that 
||f *^|| = = 1. Since, for each w'' we have e-^^w'^ = e^'=^{tv^ + w^) then 

||e-^*^i;'=|| = \e>^^^\\\{tv'' + w'')\\ > e^"^^*(t - 1). 

For any constant C > we take t > 2C + 1 and for this t we take big enough k such 
that e^'^'*'** > |. Then we have: 

||e^V|| > ^(2C + 1-1) = C 

and we conclude that He-^*]! > C for t > 2C + 1. Therefore, the family of exponents e-^* 
is not uniformly bounded and because of Banach-Steinhaus theorem there exists x G D{A) 
such that ||e'^*x|| — > oo when t +oo. 

Thus, the system flS.lOSp is unstable when s = 1. The last completes the proof of the 
proposition. □ 



6 Stabilizability analysis 

In this section we prove Theorem 11.21 for control systems (11. 8p with dety4_i = 0. It is 
convenient to reformulate this theorem as follows. 

Theorem 6.1 (on stabilizability). Let bi,...,bp G C" be the columns of the matrix B. 
Assume that the following four conditions are satisfied. 

(1) All the eigenvalues of the matrix A_i satisfy |/i| < 1. 

(2) All the eigenvalues /i G ai are simple. 

(3) Y7i=i l(^i)l/)c" 7^ for all i = 1, . . . ,p and all vectors y satisfying y G KerA^(A) for 
roots A of the equation det A^(A) = 0, such that ReA > 0. 

(4) J2^=i \ {bi,ym)c" 7^ for all vectors ym of the matrix A*^-^^, corresponding to eigenvalues 
'Prn> f^m ^ o"i cind for all i = 1, . . . ,p. 

Then there exists a regular control u = Tx of the form 

/o /-O 
F2{9)zt{9) de + J F3{9)zt{9)d9 (6.109) 

where x = iy,z{-)) G V{A), F2(-),i^3(-) G L2(-l, 0; C"^^). And this control stabilizes the 
system ITR) . i.e. V{A) = V{A + BJ^) and e^^-^+^-^^xo ^ Q as t oo for all xq G Mg. 

The controllability conditions (3)-(4) of Theorem 16.11 are equivalent to (3)-(4) of The- 
orem [T2] (for more details concerning these conditions see [12] )• We also note that regular 
stabilizability for a particular case of the control systems (11.81) had been considered in [TH] . 
Before giving the proof, let us discuss the conditions (1) and (2). Since the regular feed- 
back does not change the matrix A^i we need the assumption ct{A^i) C {/i : < 1}. 
Moreover, taking into account the results of Theorem 15.11 on strong stability, we conclude 
that by means of a regular feedback it is possible to stabilize the system in the case when 
the algebraic multiplicity of each eigenvalue ft E o-i equals to 1. In this case the closed loop 
system will be asymptotically stable if and only if all the eigenvectors of ^ + BJ-' are in 
the left half-plane. Thus, the problem of regular stabilizability for such systems consists in 
assigning the spectrum of the system in the left half-plane. The most intensional problem 
in the situation appears when an infinite number of the eigenvalues located "close" to the 
imaginary axis belongs to the right half-plane. This means, in particular, that cxi ^ 0. 
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Proof of Theorem 16. il Let us construct the decomposition of the spectrum of A introduced 
in Section [3l2] by ^M>: 

ct(^) = Ao(^)UAi(^)UA2M), (6.110) 

where the subsets are given by fl3.42p . Since Ao(^) C {A : ReA < —e}, e > then our 
aim is to construct the feedback which moves the eigenvalues of Ai(^) and A2(.4.) to the left 
half-plane. 

We begin with moving the spectral set A2(.4.). The spectral decomposition of the state 
space M2 = M2 © M2 © M| (Theorem 13.41) . corresponding to (16.1101) . allows us to rewrite 
the system (II. 9p as 

Xo = Aoxo + BqU, Xo e M° = M^{A) 

xi = Aixi + Biu, xi G = Mf{A) (6.111) 
xi = A2X2 + B2U, xi e Ml = M^{A), 

where ^0 = -^U/"? Ai = A\m^, A2 = A\]yp and Bi are the projections of B onto subspaces 

The spectrum A2(^) of the operator A2 defined on the finite-dimensional subspace 
M| belongs to the right half-plane. Due to the assumption (3) of the theorem all the 
eigenvalues of A2 are controllable and, thus, they may be assigned arbitrary to the half- 
plane {A : ReA < —s}. Namely, let us consider a feedback u = J^2X, T2 '■ M2 — > which 
acts as follows: 

^2X0 = T2X\ = 0, Xo G M2,Xi G M\, T2X2 = ^2X2, X2 G M^, 

where J-' : M| — ^ is a bounded operator. It is easy to see that such form of feedback 
allows to move only eigenvalues of A2(.4), i.e. Ao(^) U Ai(^) C a{A + BJ^2)- Thus, we 
conclude that there exists such that the spectrum of the closed-loop system A + BJ^2 is 
of the form 

(t{A + BJ^2) = fT(^)\A2(^) U A2, A2 C {A : ReA < -e}. (6.112) 

We choose the control in the form u = T2X + denote A ^=A + BTq and rewrite the 
system (II. 9p as 

x = Ax^ Bv. (6.113) 

We emphasize that due to the form of the feedback JF2 the infinite-dimensional system (I6.113p 
corresponds to the neutral type system (II. 8p with the same matrix A_i. 

Let us now construct the decomposition of the spectrum (13.411) for the operator A. Due 
to (I6.112p . we obtain A2(^) = and, thus, 

(r{A) = Ao{A) U Ai(l), Ai{A) = AiiA). 

Applying Theorem 13. 4[ we construct the spectral decomposition of the state space M2 = 
© Mj. The operator model 06.1131) may be rewritten as follows: 

= AoXo + Bov, Xo G MO = M^{A) ,^ 
-- Aixi + Biv, xi G Ml = Ml (A), 

where ^0 = -^Ia/O; Ai = A\m'^ and Bi are the projections of the vectors of B onto subspaces 
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Due to Theorem 14.41 the restriction of the resolvent R{X, A)\]y.fO is uniformly bounded on 

the set {A : ReA > 0} and, therefore, the semigroup {e^'^\M°}t>o is exponentially stable. 

To stabilize the second equation we apply the approach introduced in p2j which is based 
on the abstract theorem on infinite pole assignment. Below, for the sake of completeness and 
due to the specific form of the operator Ai, we give a simplified formulation of the mentioned 
theorem (Theorem 16. 5p . 

The theorem on infinite pole assignment holds for a single input system. However, as 
it is shown in [22], the multivariable case may be reduced to the single input case by the 
following considerations (see [22] and also the classical result on finite-dimensional control 
systems in [HO]). 

Let us consider the vector b given by 



b \ ( cihi + ... + Cphp \ ^ f B ^ 



c = Be, c 



(6.115) 



where {bi, . . . , bp} C C" are the columns of the matrix B. The exist numbers q G C such 
that the following relations hold (Lemma 16.61) : 

(a) {b, ym)c" 7^ for all eigenvectors i/m of the matrix A*_-^ corresponding to the eigenvalues 
7V, where /x^ G (Xi, m = 1, . . . , ^i. 

(b) (6, t/^)c" 7^ for all such y'^ which satisfy G KerA*(A^) and ReA^ > 0. 

Let us denote the projection of b onto the subspace by bi and consider the single 
input system 

Xi = AiXi + hiu, mgC. (6.116) 

Due to the construction, the eigenvalues of Ai are simple and each eigenvalue A^, 
m = 1, . . . , ii, \k\ > Ni belongs to the circle L^{r^^^) centered at \^ = i(arg/i.m + 27r/i;). The 
corresponding eigenvectors ip'^ form a Riesz basis (Proposition 12. 7p . 

Thus, the single input system (16.1161) satisfies the conditions (H1)-(H4) of Theorem 16. 5[ 
According to that theorem the spectrum Ai(^) of the operator Ai may be moved by the 
regular feedback (16.1091) as follows. 

Let us chose scalars A^ such that ReA^ < and which are located inside the circles 
L'^{r^''^). There exists a feedback u = T\x, T\ : M2 C such that A^ are eigenvalues 
of the operator A\ + biJ^i. Taking into account (I6.115p . we define the feedback v = J-'iX, 
Ti:M2^ CP as follows 

r-, def 

l'lX = CJ-iX. 

Due to the theorem on stability the semigroup {Q^^'^^^'^^'^\Ml]t>o is asymptotically stable. 
Thus, the feedback 

U = J^iXi + J^2X2 

transforms the original system into a system where all the conditions of Theorem 15.11 on 
asymptotic stability are verified. The last completes the proof of Theorem 16.11 □ 

Remark 6.2. We would like to emphasize that in the present paper to prove the result on 
stabilizability for the case detyl„i = we have contributed the ideas which are mainly of 
technical character: the direct decomposition of the state space and the proof of the resolvent 
boundedness on some subspace. However, the main contribution from the stabilizability point 
of view is the abstract theorem on infinite pole assignment which had been proved in JM 
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Remark 6.3. The stabilizability of the restriction of the system onto the subspace M|(^) 
follows also from some classical results. Let us denote by Ts a rectifiable, simple, closed curve 
which surround the spectral set A2{A) and C {A : ReA > e}. By P2 = f R{A, A) dA 

we denote the spectral projector. The subspaces P2M2 and (/ — P2)M2 are A-invariant 
and the spectrum of the restriction (/ — P2)^2 belongs to the half-plane {A : ReA < e}. 
According to JW^ the spectral set h.2{A) may be assigned arbitrary to the half-plane {A : 
ReA < e} by means of a finite rank input operator (under some controllability conditions) . 
The development of this approach was given in fL^ and (see e.g. as a survey). 

Remark 6.4. Let us discuss the assumptions (3)-(4) of Theorem \l.^ - 

(3) rank(A(A), B) = n for all A : ReA > 0. 

(4) rank(/i/ — A^i, B) = n for all fi E ai. 

The assumptions (3) may include an infinite number of relations (obviously, it should 
be verified only for A which are eigenvalues of A, i.e. det A_4(A) = 0; there may exist an 
infinite number of eigenvalues of A belonging to the right half-plane). However, only a finite 
number of these relations have to be verified. More precisely, there exists M > such that 
for any eigenvalue A such that |ImA| > M the condition (3) follows from the condition (4). 

Proof. Let A^i be of the form fl2.16p and we prove the proposition for eigenvalues located 
inside the circles L^(r*^^)), i.e. for A^ (for all other indices m = 1, . . . , £1 the idea of the proof 
remains the same). The matrix A^i — fiil = A^i — e~^'^I is of the form: 

A_i-^i/=(^° lGd"-^)^("-i), detl^O. 

Due to the assumption (4) there exists i G l,p such that the column b = {bu, . . . ,bni)'^ of 
thee matrix B possesses a nonzero first component: bu 7^ 0. Besides, we note that since 
det A ^ then there exists the unique a = (02, . . . , an)"^ £ C"~^ such that b = Aa, where 
b = {b2i,...,bnif eC"-\ 

Let us analyze the structure of A(Aj^), we rewrite it as follows: 

A(At) = Ate-^^(A_i- /ill + 

where 

D, = (e^^ _ e^^)/ + e^^ e^^^ (^A2{9) + j^^A,{9)^ d9. 

From the last formula we conclude, that 0, when k 00. 

Let us show that there exists G N such that for any k : |fc| > iV: 

rank(A(At), b) = n, (6.117) 

which is equivalent to the statement of the proposition. 

Since det A(A^) = then det{A_i- fij + Dk) = 0. Let us denote by Dk = {(A)y}"j=2- 
Since det A ^ 0, then there exists N E N such that for any k : > we have: det (A + 
Dk) 7^ 0. Therefore, there exists the unique vector ak G C" such that b = {A-\- Dk)ak. From 
the last we conclude 

= {A + Dk)-'b= {I + A''Dk)-'A-'b= {I + A-'Dkr'a. 

Since (/ + A^^D^)^^ I when — > 00, we conclude that — > a. 
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If we suppose that fl6.117l) does not hold, then 61 = bn allows the following representa- 
tion: 

hi = {ak)2{Dk)l2 + {0ik)n{Dk)ln- 

However, since (-Dfc)ii 0, then the right-hand side of the last relation tends to zero when 
k 00. We have come to the contrary, which completes the proof of the proposition. □ 

As it was mentioned above, for the sake of completeness we give further the formulations 
of the results on infinite pole assignment from [22]. We formulate them in a form, which 
takes into account the specific form of the operator Ai. 

Theorem 6.5 (On infinite pole assignment [22]). Let H be a complex Hilbert space, A 
be an infinitesimal generator of a Co-semigroup in H , and the control system if given by 
X = Ax + Bu, X G 'D{A) C H. Let fii, . . . , fig be some nonzero complex numbers and we 
introduce complex numbers 

= In l/i^l + i(arg/im + ^nk), m = 1, . . . , £, A; G Z, 

and the circles L^(r*^'^)) centered at with radii r^'^^ satisfying the relation ^ (r*^^))^ < 00. 
Let the following assumptions hold: 

(HI) The spectrum of A consists only of eigenvalues which are located in the circles L^{r^^'>). 
Moreover, all the eigenvalues of A are simple (i.e. its algebraic multiplicity equals 1) 
and there exists A?"! G N such that for any k : \k\ > Ni the total multiplicity of the 
eigenvalues contained in the circles L^{r^''^) equals 1. 

(H2) The corresponding eigenvectors, which we denote by ip'^, constitute a Riesz basis in H . 

(H3) The system x = Ax + Bu is of a single input, i.e. the operator B : C H is the 
operator of multiplication by h E H. 

(H4) We assume the following controllability condition: b is not orthogonal to eigenvectors 
ip^ of the operator A* : (h^ip^ 7^ and 

lim A;|(b,7/;^)| = c^, < c„ < +00. (6.118) 

fc— >co 

Then there exists N2 > Ni such that for any family of complex numbers Aj^ G L^{r^^'>), 
\k\ > N2 there exists a linear control T : T){^A) C, such that 

(1) the complex numbers A^ are eigenvalues of the operator A + BT ; 

(2) the operator BJ-' : T^IA) ^ H is relatively A-bounded. 

Lemma 6.6 ([22j). If the assumptions (3) and (4) of Theorem \6. 1\ holds, then there exists 
a vector b G ImB, say b = Ci6i + . . . + Cpbp, Ci G C, such that the following relations hold: 

(a) {b,ym)c" 7^ for all eigenvectors ym of the matrix A*_^ corresponding to the eigenvalues 
/i;^, where /i^ G cti, m = 1, . . . , ^1. 

(b) (6, 2/m)c" 7^ for all such which satisfy G KerA*(A^) and ReA^ > 0. 
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Proof. Let us prove the relations (a). Consider the subspaces 

Lm = {y- (y,ym)c" = 0}, m=l,...,£i. 

These subspaces are of dimension n — 1. Let us consider also = H ImB. By the 
assumption (4) we have: dim Mm < p. Indeed, if dimM^ = p = dim ImB, then B*y = 
and the condition (4) is not satisfied. Thus, each Mm is nowhere dense in Imi? and due to 
Baire theorem we have: 

ii 

[j ImB. 

m=l 

The last means that there exists h G Imi? such that h ^ Lm for all m = 1, . . . i.e. the 
relations {b,ym)c" 7^ hold. 

The inequalities (b) follows from the relation 

\{b,yt)c^\^CU\k\ + ir' 
which is proved in [22j independently of the condition det^l-i 7^ 0. □ 

7 Conclusion and perspectives 

In the present paper we have generalized the results on strong asymptotic non-exponential 
stability and regular stabilizability for the case of mixed retarded-neutral type systems. 
The proofs of the mentioned generalizations are technically complicated and requires subtle 
estimates. We combine the Riesz basis technic with the analysis of the boundedness of the 
resolvent on some ^-invariant subspaces. 

As a perspective, we consider systems with the difference operator K given by 

r 

Besides, the dilemma of item (iii) of Theorem 11.11 on stability may be investigated more 
precisely. 
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